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Abstract. 

We analyze the asymptotics of the Wigner 3,j'-symbol as a matrix element 
connecting eigenfunctions of a pair of integrable systems, obtained by lifting 
the problem of the addition of angular momenta into the space of Schwingcr's 
oscillators. A novel element is the appearance of compact Lagrangian manifolds 
that are not tori, due to the fact that the observables defining the quantum 
states are noncommuting. These manifolds can be quantized by generalized Bohr- 
Sommerfeld rules and yield all the correct quantum numbers. The geometry of the 
classical angular momentum vectors emerges in a clear manner. Efficient methods 
for computing amplitude determinants in terms of Poisson brackets are developed 
and illustrated. 



PACS numbers: 03.65.Sq, 02.20.Qs, 02.30.1k, 02.40.Yy 

1. Introduction 

This article is a study of the asymptotics of the Wigner 3j-symbol from the standpoint 
of semiclassical mechanics, that is, essentially multidimensional WKB theory for 
integrable systems. The principal result itself, the leading asymptotic expression for 
the 3j-symbol, has been known since Ponzano and Regge (1968). Nevertheless our 
analysis presents several novel features. One is the exploration of Lagrangian manifolds 
in phase space that are not tori (the usual case for eigenstates of integrable systems) . 
Instead, one of the states entering into the 3j-symbol is supported semiclassically on 
a Lagrangian manifold that is a nontrivial 3-torus bundle over 5*0(3). This manifold 
can be quantized by generalized Bohr-Sommerfeld rules, whereupon it yields the exact 
eigenvalues required by the quantum 3j-symbol, as well as the correct amplitude and 
phase of its asymptotic form. This unusual Lagrangian manifold arises because the 
quantum state in question is an eigenstate of a set of noncommuting operators. Other 
novel features include the expression of the asymptotic phase of the 3j-symbol in terms 
of the phases of Schwinger's harmonic oscillators and the determination of stationary 
phase points by geometrically transparent operations on angular momentum vectors 
in three-dimensional space. Yet another is the representation of multidimensional 
amplitude determinants as matrices of Poisson brackets. Representations of this type 
have been known for some time, but they are generalized here to the case of sets of 
noncommuting operators. The final result is a one-line derivation of the amplitude of 
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the asymptotic form of the 3j-symbol. Similarly brief derivations are possible for the 
amplitudes of the 6j- and 9j-symbols. 

In addition our analysis of the 3j-symbol may prove to be useful for the asymptotic 
study of the 3nj-symbols for higher n. The leading order asymptotics of the 6j-symbol 
were derived by Ponzano and Regge (1968), but the understanding of the asymptotics 
of the 9j-symbol is still incomplete. These symbols are important in many applications 
in atomic, molecular and nuclear physics, for example, the 9j-symbols are needed in 
atomic physics to convert from an LS'-coupled basis to a jj-coupled basis. These 
symbols are all examples of closed spin networks, of which more elaborate examples 
occur in applications, each of which presents a challenge to asymptotic analysis. 
Moreover in recent years new interest in this subject has arisen from researches into 
quantum computing (Marzuoli and Rasetti, 2005) and quantum gravity, where new 
derivations of the asymptotics of the Wigner 6j-symbol have been produced as well 
as generalizations to other groups such as the Lorentz group. The 3nj-symbols and 
their asymptotics have also been used recently in algorithms for molecular quantum 
mechanics (De Fazio et al 2003 and Anderson and Aquilanti, 2006), which exploit the 
connections with the theory of discrete orthogonal polynomials (Aquilanti et al 1995, 
2001a, b and references therein). 

The asymptotic formula for the 3j-symbol is closely related to that for the 6j- 
symbol, being a limiting case of the latter. These were first derived by Ponzano 
and Regge (1968), using intuitive methods and building on Wigner's earlier result for 
the amplitude of the 6j-symbol (Wigner, 1959). Later Neville (1971) analyzed the 
asymptotics of the 3j- and 6j- symbols by a discrete version of WKB theory, applied 
to the recursion relations satisfied by those symbols, without apparently knowing of 
the work of Ponzano and Regge. His formulas are not presented in a particularly 
transparent or geometrical manner, but appear to reproduce some of the results of 
Ponzano and Regge. The formula for the 3j-symbol (in the form of Clebsch-Gordan 
coefficients) was later derived again by Miller (1974), who presented it as an example 
of his general theory of semiclassical matrix elements of integrable systems. Miller 
called on the fact that the phase of the semiclassical matrix element is a generating 
function of a canonical transformation, and used the classical transformation that most 
obviously corresponds to the quantum addition of angular momenta to reconstruct the 
generating function. The method leads to a difficult integral, which, once done, yields 
the five terms in the phase of the asymptotic formula for the 3j-symbol. Somewhat 
later Schultcn and Gordon (1975a, b) presented a rigorous derivation of the Ponzano 
and Regge results for the 3j- and 6j-symbols, using methods similar to those of 
Neville but carrying them out in a more thorough and elegant manner. Schulten and 
Gordon also provided uniform approximations for the transition from the classical to 
nonclassical regimes, work that has recently been reanalyzed (Geronimo et al , 2004) 
and extended to non-Euclidean and quantum groups (Taylor and Woodward, 2005). 
Somewhat later Biedenharn and Louck (1981b) presented a review and commentary 
of the results of Ponzano and Regge, as well as a proof based on showing that the 
result satisfies asymptotically a set of defining relations for the 6j-symbol. More 
recently the asymptotics of the 3j-symbol was derived again by Reinsch and Morchcad 
(1999), working with an integral representation constructed out of Wigner's single- 
index sum for the Clebsch-Gordan coefficients. About the same time, Roberts (1999) 
derived the Ponzano and Regge results for the 6j-symbol, using methods of geometric 
quantization. Finally, Freidel and Louaprc (2003) presented a derivation of the 
asymptotic expression for the square of the 6j-symbol, based on an analysis of an 



Semiclassical analysis of Wigner 3j -symbol 



3 



SU (2) path integral. This work was part of a larger study of generalizations of the 6j- 
symbol to other groups (for example, the lOj-symbol) that are important in quantum 
gravity. See also Barrett and Steele (2003) and Baez, Christensen and Egan (2002). 

There are many variations on the calculation of the asymptotic forms of the 3nj- 
symbols that have been considered by different authors. There are asymptotic forms 
inside and outside the classically allowed regions; uniform approximations connecting 
two or more of these regions; asymptotic forms when only some of the quantum 
numbers are large and others small; and higher order terms. Ponzano and Regge 
(1968) covered many of these issues, while Reinsch and Morehead computed some 
higher order terms. 

The outline of this paper is as follows. In Sec. [2] we review the semiclassical 
mechanics of integrable systems in the generic case that one has sets of commuting 
observables, drawing attention to an expression for the amplitude determinant in terms 
of Poisson brackets. In Sec. [3]we review the Schwinger model for representing angular 
momentum operators in terms of harmonic oscillators. This model allows us to express 
angular momentum eigenstates in terms of wave functions on R™, which we use in Sec. [4] 
to express the 3j-symbols in terms of scalar products of such functions. In Sec. Owe 
study the Schwinger model from a classical standpoint, in which an important element 
is the reduction of the Schwinger phase space (the "large phase space" ) by the torus 
group T 3 , producing the Poisson manifold R 3 x I 3 x R 3 ( "angular momentum space" ) 
and the reduced phase space S 2 x S 2 x S 2 (the "small phase space"). In Sees. [6] and [7] 
we study the two Lagrangian manifolds that support the states whose scalar product is 
the 3j-symbol. One is a conventional invariant torus (the "Jm-torus"), but the other, 
what we call the "Wigner manifold," is compact and Lagrangian but not a torus. This 
manifold supports Wigner's state of zero total angular momentum that enters into the 
definition of the 3j-symbols. In Sees. [8] we study the intersections of the jm-torus and 
the Wigner manifold, which are the stationary phase points of the 3j'-symbol, and show 
how these can be found by elementary geometrical considerations in three-dimensional 
space (that is, by rotating angular momentum vectors). The intersection of the two 
manifolds turns out to be a pair of 4-tori. In Sec. [5] we compute the action integrals 
along the respective Lagrangian manifolds to points on the two 4-tori, whose difference 
is the Ponzano and Regge phase of the 3j-symbol. In Sec. [10] we apply generalized 
Bohr-Sommcrfeld quantization to the jm-torus and the Wigner manifold, a standard 
procedure for the jm-torus, although it leads in an interesting way to the extra 1/2 in 
the classical values representing the lengths of the angular momentum vectors. This 
extra 1/2 was guessed by Ponzano and Regge and Miller and derived systematically 
by Schulten and Gordon, Reinsch and Morehead and by us, although it is missing 
from the results of Roberts. In our work it is essentially a Maslov index. In Sec.fTTIwe 
generalize known expressions for the amplitude determinant of semiclassical matrix 
elements of integrable systems in terms of Poisson brackets to the case of collections of 
noncommuting observables (whose level sets nevertheless are Lagrangian). The result 
allows us to compute the amplitude of the 3j-symbol as a 2 x 2 matrix of Poisson 
brackets. We then put all the pieces together to obtain the final asymptotic form. 
Finally, in Sec.[12]we present some comments on the work, prospects for further work, 
and conclusions. 
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2. Semiclassical wave functions for integrable systems 

The semiclassical mechanics of integrable systems is well understood (Einstein, 1917; 
Brillouin, 1926; Keller, 1958; Percival, 1973; Berry and Tabor, 1976; Gutzwiller, 1990; 
Brack and Bhaduri, 1997; Cargo et al , 2005a, 2005b). Here we summarize the basic 
facts, some of which require modification for our application. 

We consider the quantum mechanics of a particle moving in W 1 (with wave 
function ip(x\, . . . , x n ) and Hilbert space L 2 (M. n )). We speak of an integrable system if 
we have a complete set of commuting observables {Ai, . . . , A n } acting on this Hilbert 
space. We use hats to distinguish quantum operators from classical quantities with a 
similar meaning. Sometimes the Hamiltonian is one of these operators or a function 
of them, but in our application there is no Hamiltonian, or, rather all the A^s are 
Hamiltonians on an equal footing. These operators may be converted into their 
classical counterparts by the Weyl transform (Weyl 1927, Wigner 1932, Groenewold 
1946, Moyal 1949, Voros 1977, Berry 1977, Balazs and Jennings 1984, Hillery et al 
1984, Littlejohn 1986, McDonald 1988, Estrada et al 1989, Gracia-Bondfa and Varilly 
1995 and Ozorio de Almeida 1998). The Weyl transforms (or Weyl "symbols") of 
these operators are functions on the classical phase space R 2 ™, that is, functions of 
(xi, . . . , x n ;pi, . . . ,p n )- They are normally even power series in h, as we assume, of 
which the leading term is the "principal symbol." We denote the principal symbols of 
{Ai, . . . , A n } by {Ai, . . . , A n } (without the hats). In view of the Moyal star product 
representation (Moyal 1949) of the vanishing commutators [Ai, Aj] — 0, the principal 
symbols Poisson commute, {Ai,Aj} = 0, thus defining a classically integrable system 
(Arnold 1989, Cushman and Bates 1997). (We use curly brackets { } both to denote a 
set and for Poisson brackets.) Then according to the Liouville- Arnold theorem (Arnold 
1989), the compact level sets of {A±, . . . , A n } are generically n-tori. The Hamiltonian 
vector fields generated by the Ai are commuting and linearly independent on the tori; 
thus the tori are not only the level sets of the Ai, they are also the orbits of the 
Abclian group generated by the corresponding Hamiltonian flows. One can define an 
action function S on a torus as the integral of ^2nP% dxi relative to some initial point; 
it is multivalued because of the topologically distinct paths going from the initial to 
the final point, but otherwise is independent of the path. 

Let Ai = ai be one of these tori (Ai are the functions, a* the values). The torus 
has a projection onto configuration space defining a classically allowed region in that 
space; the inverse projection is multivalued. The function S may be projected onto 
configuration space, defining a function we shall denote by Sk(x, a) (where for brevity 
x and a stand for (xi, . . . , x n ) and (ai, . . . ,a n ), etc). Here k labels the branches 
of the inverse projection; function Sk has an additional multivaluedness due to the 
choice of contour connecting initial and final points on the torus. Then as explained 
by Arnold (1989), S(x, a) is the generating function of the canonical transformation 
(x,p) — > (a, A), where a = (ai, . . . , a n ) is the set of angle variable conjugate to the 
conserved quantities (Ai, . . . , A n ). Action variables may be defined in the usual way as 
(l/27r) jpdx around the independent basis contours on the torus; these are functions 
of the Ai and their conjugate variables are the angles that cover the torus once when 
varying between and 27r. Sometimes however it is more convenient to work with the 
Ai instead of the actions (the Ai are not necessarily actions, and their flows are not 
necessarily periodic on the torus). 

Tori arc quantized, that is, associated with a consistent solution of the 
simultaneous Hamiltonian- Jacobi and amplitude transport equations for the operators 
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Aj, only if they satisfy the Bohr-Sommerfeld or EBK quantization conditions, 
discussed in Sec. [TO] Associated with a quantized torus is a semiclassical wave function 
in configuration space, which in the classically allowed region is given by 

1&(a0 = (x\a) = ^=^2\n k \ l / 2 eMi[Sk(x,a)- f i k 7r/2}}. (1) 

The meaning of this formula is the following. First, here and below we set h = 1. 
Next, given the point x in the classically allowed region, its inverse projection onto 
the quantized torus is a set of points indexed by k. We assume the projection is 
nonsingular at these points (we are not at a caustic). The phase S k (x,a) is the 
integral of p dx from a given initial point on the torus to the fc-th point of the inverse 
projection, and [i k is the Maslov index (Maslov 1981, Mishchenko et al 1990, de 
Gosson 1997) of the same path. The amplitude determinant Q k is given by 

n k = det = [det{x i} A,}}- 1 , (2) 

where in the second form the Poisson brackets are evaluated on the fc-th branch of 
the inverse mapping from x to the Lagrangian manifold. The amplitude determinant 
is a density on configuration space (to within the semiclassical approximation, the 
probability density corresponding to a single branch), which is the projection onto 
configuration space of a density on the torus. The latter density is required to be 
invariant under the Hamiltonian flows generated by the Aj (this is the meaning of 
the amplitude transport equations for the Ai); in terms of the variables a, conjugate 
to the Ai this means that the density is constant (it is the n-form dct\ A ... A da n ). 
Finally, the quantity V in ^ is the volume of the torus, measured with respect to 
this density. If the Ai are action variables, then V — (27r) n . The overall phase of the 
wave function (its phase convention) is determined by the choice of the initial point 
on the torus. 

Now let {Ai, . . . , A n } and {B\, . . . , B n } be two complete sets of commuting 
observables, with principal symbols Ai and Bi, conjugate angles etj and (3i and action 
functions Sa{x,o) and Ss{x,b), and let a and b refer to two quantized tori (an A- 
torus and a -B-torus). We assume initially that the two sets of functions Ai and Bi are 
independent. We compute (b\a) as an integral of the wave functions over x, evaluated 
by the stationary phase approximation. The stationary phase points are geometrically 
the intersections of the A-torus with the -B-torus. Generically the two tori intersect 
in finite set of isolated points that we index by k, denoting the corresponding a and j3 
values by a k and (3 k . (For given k, a k and k refer to the same point in phase space.) 
Then the result is 

(b\a) = i-=i=^|^7 fe | 1 /2exp{^[5 A K)-5B(/3 fc )-M^/2]}. (3) 
VVaVb 

Here Va and Vg are the volumes of the respective tori, as in (Q}, and the actions Sa 
and Sb are considered functions of the a or (3 coordinates on the respective tori. 

As shown by Littlejohn (1990), the amplitude determinant Q k can be written in 
terms of the Poisson brackets of the observables Ai and Bi, 

n k = [det{A i ,B j }}~ 1 . (4) 

The Maslov index /i k in is not the same as in (JXJ) . 

Another case considered by Littlejohn (1990) is the one in which some of the 
Ai are functionally dependent on some of the Bi. For this case it is convenient 
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to assume that the first r of the two sets of variables A and B are functionally 
independent, while the last n — r are identical, so that A = {Ai, . . . , A r ,A r+ i, . . . , A n } 
and B = {B±, . . . , B r , A r+ i, . . . , A n }. Then the stationary phase points are still the 
intersections of the two n-tori, but now the intersections are generically a finite set of 
isolated (n — r)-tori, upon which linearly independent vector fields are the Hamiltonian 
vector fields associated with the (A r +i, . . . , A n ). Such an (n — r)-torus is the orbit of 
the Abelian group action generated by the corresponding Hamiltonian flows. In this 
case we find 



where now Vk is the volume of the /c-th intersection (an (n — r)-torus on which the 
volume measure is da r +i A ... A da n ), and where now the amplitude determinant f2fc is 
still given by ([4]), except that it is understood that only the first r of the A's and £Ts 
enter (thus, it is an r x r determinant instead of an n x n one). The phase difference 
Sa — Sb for branch k can be evaluated at any point on the (n — r)-torus which is 
the intersection, since the integral of pdx back and forth along a path lying in the 
intersection vanishes. 

3. The Schwinger model 

The Schwinger (or SU(2) or boson) model for angular momenta is explained well in 
Schwinger's original paper (reprinted by Biedenharn and van Dam (1965), the original 
1952 paper being unpublished), and reworked in an interesting way by Bargmann 
(1962). For further perspective see Biedenharn and Louck (1981a) and Smorodinskii 
and Shelepin (1972). Introductions are given by Sakurai (1994) and Schulman (1981). 
Here we define the notation for the Schwinger model and emphasize some aspects that 
will be important for our application. 

In the Schwinger model each independent angular momentum vector is associated 
with two harmonic oscillators. We shall refer to the lj-, 3j-, etc models, depending on 
how many independent angular momenta there are. The number of j's in the model 
is not necessarily the number of j's in the Wigner symbol; for example, Miller (1974) 
used a 2j-model to study the Clebsch-Gordan coefficients, essentially the 3j-symbols. 

We start with the lj-model, for which there are two harmonic oscillators indexed 
by Greek indices /x, u, . . . = 1, 2. (These are just labels of the two oscillators; sometimes 
other labels such as 1/2, —1/2 are more suitable.) The wave functions are ip(xi,X2) 
and the Hilbert space is H = L 2 (R 2 ). We write = (1/2) (i^ + p2) for the two 

oscillator Hamiltonians, and we define H = ^ M H^. The eigenvalues of H are n + 1, 
with n = 0, 1, . . ., and energy level E n is (n + l)-fold degenerate. We introduce usual 
annihilation and creation operators a M = (x^ + ip M )/v / 2, = (x^ — ipf^/V^, omitting 
the hats on the a's and a^'s since these will always be understood to be operators. We 
define operators 





(6) 



and 




(7) 
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where a 1 is the i-th Pauli matrix. Here and below we use indices ■ = 1,2,3 

(or x,y,z if that is more clear) to denote the Cartesian components of a 3-vector. 
Notice that I and Jj are quadratic functions of the x's and p's of the system. The 
eigenvalues of / are n/2 for n = 0, 1, These operators satisfy the commutation 

relations [I,Ji] = and [Jj, Jj] = i^2 k £ijkJk- We also define J = J2i^h so that 
[I, J 2 ] = and [Ji, J 2 ] = 0. It avoids some confusion with indices to always denote 
the square of a vector by a bold face symbol, as we have done here. We note the 
important operator identity J 2 = /(/ + 1), expressing the quartic operator J 2 as a 
function of the quadratic operator /. 

From this identity and the known eigenvalues of I it follows that the eigenvalues 
of J 2 are (n/2) [(n/2) + 1], for n = 0,1,..., which leads us to identify n/2 with 
j = 0, 1/2, 1, . . ., the usual angular momentum quantum number. The n-th (or j-th) 
eigenspace of H or / is (2j + l)-dimensional, and so must contain a single copy of 
the j-th irrep of SU(2). Each irrep (both integer and half-integer values of j) occurs 
precisely once in the Hilbert space TL. We denote these subspaces by Hj, and write 
H = ^ • (BHj. The standard basis in Hj is the eigenbasis of J z = \{H\ — H 2 ), with 
the usual quantum number m, so that if n M are the usual quantum numbers of the 
oscillators H^, then n\ = j + to, n 2 = j — m. The simultaneous eigenstates of J 2 and 
J z are \jm) or |nin 2 ). 

In the TVj-model we index the angular momenta with indices r, s, . . . = 1, . . . , N. 
The oscillators are now labelled H r ^ with coordinates and momenta x r ^ and ff r ^ 
and annihilation and creation operators a rfl and aj^. The wave functions are now 
tp(xn,Xi2,X2i, ■ ■ ■ ,x N2 ) and the Hilbert space is L 2 (R 2N ). We define operators 





i v t 


(8) 


Jri 


— 2 Z-^i a " r t J - a t* v a ' ri/1 


(9) 


T 2 


- Vi 2 

i 


(10) 


J t 


= ^Jri, or J = ^ J r , 

r r 


(11) 


j 2 




(12) 



i 



most of which are obvious generalizations from the case N = 1. These satisfy the 
identity 

J 2 =/ r (ir + l), (13) 

and the commutation relations [I r , J s i\ — [I r , J 2 ] = 0. Each angular momentum vector 
J r also obeys the standard commutation relations among its components and square, 
which we omit, as does any sum of these angular momenta (partial or total). 

The angular momenta generate an action of [SU(2)] N on the Hilbert space (one 
copy for each pair of oscillators). Here we discuss only the simultaneous rotation of 
all oscillator degrees of freedom by the same element of 5/7(2), which is generated by 
the total angular momentum, but partial rotation operators can also be defined and 
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are useful. We begin with the commutation relations, 

[Ji,a r ^]=--'Y^a % llv aru, (14) 

V 

which define the transformation properties of the operators a r(tt , a\ „ under infinitesimal 
rotations. We define a finite rotation operator in axis- angle or Euler angle form by 

U(n,6) =exp(-i0n-J), (16) 

U(a, 0, 7 ) = U(z, a)U(y, (3)U(z, 7), (17) 

where n is a unit vector defining an axis and 9 an angle of rotation about that axis, 
and where x, y and z are respectively the unit vectors along the three coordinate axes. 
The U operators form a faithful representation of SU(2). 

We use the symbol u(n, 9) or u(a, f3, 7) for the 2x2 matrices belonging to SU (2), 
in axis-angle or Euler angle parameterization (not to be confused with the U operators 
that act on the Hilbert space of the 2 N oscillators) . Thus 

u(n,6) = exp(-i6>n- a/2) = cos6»/2 - in- cr sin 0/2. (18) 

The exponentiated versions of equations (fT4"|) and (fT5]) are 

O r ^ L U ^ ^ tL^iv ^rvi 



(19) 

rfi / j 



where both U and u have the same parameterization. In the language of irreducible 
tensor operators the pair of operators (a r i, a r 2) transforms as a spin-1/2 operator. 

Similarly, vector operators are the angular momenta themselves, which satisfy 
the conjugation relations, 

tfJnU = ^2RijJ rj , (20) 

j 

where R is the 3x3 orthogonal rotation matrix with the same axis and angle as U. 
The relation between R and u (with the same axis and angle) is 

R4o=\tx{U^a i Ua j ). (21) 

This is the usual projection from SU(2) to 5*0(3), in which the inverse image of a 
given R € 5*0(3) is a pair (u, —u) in SU(2). 



4. The Wigner 3j-symbols in the Schwinger Model 

We now define the 3j-symbols in the context of the Schwinger model. We take 
the 3j-model, N — 3. One complete set of commuting observables on the 
Hilbert space H. <g> 7i <g> H is (Ii, I2, 13, Jiz, Jiz, J^z), with corresponding eigenstates 
|jij2j3Wim2m3) = |j2m2)|j3m3). Another complete set arises in the usual 

problem of addition of three angular momenta, in which we consider the values of 
j and m (the quantum numbers of J 2 and J z ) that occur in the product space 
Hj 1 <g> Hj 2 Hj 3 for fixed values of (ji, J2, J3), a subspace of H <g> H (E> H. The set 
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of the five commuting operators 1%, I3, J3, J 2 ) that arises in this way is however 
not complete (the simultaneous eigenstates in general possess degeneracies), so to 
resolve these we introduce a sixth commuting operator, conventionally taken to be 
J?2 = (Ji + J2) 2 with quantum number ji 2 (J23 or J 2 3 will also work). 

The Wigner 3j-symbols only involve the case 3 ' = 0, but we mention the others 
anyway because the foliation of the classical phase space into Lagrangian manifolds 
involves the other values. The usual rules for the addition of angular momenta show 
that if (ii,i2,j3) satisfy the triangle inequality, then there exists precisely a one- 
dimensional subspace of Hj 1 <8> Hj 2 ® Hj 3 with j = 0; if they do not, then no such 
subspace exists. If we enlarge our point of view to the full Hilbert space then 
there is an infinite dimensional subspace with 3 = 0, a basis in which is specified by 
all triplets (ji, J2, J3) that satisfy the triangle inequalities. If j = 0, then the quantum 
number m is superfluous, since m = 0; the quantum number j'12 is superfluous as well, 
since j 12 = j 3 . 

We note that if (ip\3 2 \ip) = for any state \ip), then Ji\ip) = for i = 1,2,3. 
Although the components of J do not commute and so do not possess simultaneous 
eigenstates in general, the case of a state with j = is an exception, since it is 
a simultaneous eigenstate of all three components of J with eigenvalues 0. With 
this in mind we denote the basis of states in the subspace of the full Hilbert space 
with j '■ — by \j1faj3O), where the zero vector indicates the vanishing eigenvalues 
of J. These basis states are also eigenstates of the operators J?-, for example, 

Jf2li1.72.73O) = 33(33 + WihjsO)- 

When the phase of the state \ j1j2j3O} is chosen to agree with Wigner's convention 
for the phases of the 3j-symbols, we have 



3i 32 33 

mi 7TJ2 TO3 



(ji32j3m 1 m 2 m3\j 1 j2j 3 0}. (22) 



In this manner we have expressed the 3j-symbol as a matrix element connecting 
the eigenstates of two sets of observables, (Ii, I2, is, Ji z , J2z, J3z) on the left and 
(Ii, I 2 , 13, J Xl J y , J z ) on the right. Since the second set is noncommuting, we will 
require a generalization of §5§ to compute the semiclassical approximation to the 3j- 
symbols. 



5. Classical mechanics of the Schwinger model 

The classical mechanics of the Schwinger model must be well understood in order 
to carry out a semiclassical analysis. A general reference on the classical mechanics 
of integrable systems from the modern point of view is Cushman and Bates (1997), 
where harmonic oscillators in particular are treated. 



5.1. The lj -model 

We start with the lj-model, defining two classical oscillators = (l/2)(x 2 +pfj, 
and H — ^^H^, as in the quantum case. The classical configuration space is 
M 2 and the phase space is K 4 . We introduce complex coordinates on phase space 
z ti = ( x n + iPn)/V2 an d = (x^ — ip^)/^/2, where we use an overbar for complex 
conjugation. These are the Weyl symbols of the operators a p , at. The complex 
coordinates z^, allow us to identify the phase space R 4 with C 2 , that is, knowledge 
of Z\ and z 2 allows us to find all four real coordinates (x\,X2,p\,p2), since the z's are 
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complex conjugates of the z's. As we shall see, coordinates (zi,z 2 ), arranged as a 2- 
component column vector, transform as a spinor under certain SU (2) transformations. 
Variables z M and iz^ are canonically conjugate (g's and p's respectively), so that the 
Poisson bracket of two functions / and g on phase space can be written, 

r, } = f df dg df dg 
U,3i ^\dxadpa dp^dx^ 



E( df dg df dg \ 



\dz,j, d(izn) diiz^dz^ 
The basic building blocks of the classical Schwinger model are the function 



I" M 



and the three functions 

Ji = g E! ^ ^ ( 25 ) 

for i = 1,2, 3, which define a classical angular momentum vector. We also define 
J 2 = J2i J?- These functions satisfy the identity J 2 = I 2 and the Poisson bracket 
relations {/, J;} = 0, {Ji, J,} = J2k e ijkJk, {Ji, J 2 } = 0. 

There are two groups of interest that act on the phase space R 4 or C 2 . The first 
is U(l), generated by /. Hamilton's equations for / are 

dz„ _ dl _ i _ dl _ l 

2 <ty &*„ 2 1 j 

where tp is the parameter of the orbits. These have the solution 

z„,W) = exp(-#/2) z M (0), z^ip) = exp(#/2) f M (0). (27) 

Under the /-flow, the two-component spinor [z\, z 2 ) just gets multiplied by an overall 
phase exp(—iip/2). Except for the special initial condition (z 1: z 2 ) = (0,0) (the origin 
of phase space R 4 or C 2 ), the orbits are circles with period Aw with respect to the 
variable ip. Henceforth when citing equations such as (f2"r?)) or (j2~7| we shall omit the 
second half, when it is simply the complex conjugate of the first half. 

We denote a value of / by j > 0. This is convenient notation, but in this classical 
context j is a continuous variable not to be identified with the quantum number of any 
operator (see Sec fTUj) . Except for the origin j = 0, the level set I = j (or equivalently, 
J 2 = j 2 ) is the sphere S 3 , which is foliated into circles by the action (|2"T)) . This 
foliation is precisely the Hopf fibration (Frankel 1997, Nakahara 2003), yielding the 
quotient space S 2 — S 3 /S 1 . 

The second group acting on phase space is 577(2), whose action is generated by 
the Ji. Explicitly, if n is a unit vector and an angle, then the solutions of Hamilton's 
equations 

dz u 9(n-J) ir^, , . . 

i = ^ = -2E(-^^ ( 28 ) 

and its complex conjugate are 

^(0) =^(11,0)^(0) (29) 
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and its complex conjugate. These are the obvious classical analogs of equations (O; 
notice that the period in 9 is Att. It is because of this SU(2) action that we say that 
coordinates {z\, Z2) form a spinor. This classical action of SU(2) can be understood as 
a subgroup of the classical group of linear canonical transformations, Sp{4) (Littlejohn 
1986); in general, Sp(2N) possesses a subgroup Sp(2N) n 0(2N) that is isomorphic 
to U(N), which contains the subgroup SU(N) (in this case, N = 2). When the 
symplectic matrices lying in the SU(2) subgroup are expressed in the complex basis 
(ZfijiZfj), they block diagonalize with u multiplying the z's and u multiplying the z's. 

Equation (j2"5)) defines a map (a projection) tt : M 4 (or C 2 ) — > R 3 , where M 3 is 
"angular momentum space," the space with coordinates ( J±, J2, J3). Here and below 
we use tt to denote this map or its generalization to the iVj-model. The map tt maps 
a larger space onto a smaller one, and so is not one-to-one. The inverse image of a 
point J of angular momentum space is a set of spinors that differ by an overall phase. 
It is easy to see that the definition P5|) does not depend on the overall phase of the 
spinor. Thus, the inverse image is a circle, except in the case J — when it is a single 
point (the origin of phase space C 2 or M 4 ). 

These circles are precisely the orbits of the /-flow (|27j> . Any function / that is 
constant on these circles projects onto a well defined function on angular momentum 
space. But such functions are those that Poisson commute with I, {/,J}=0. This 
includes I itself as well as the three J*. We can write such a function as /(zi, z 2 , z\, z%) 
or /(J). Now if / and g are any two such functions, then so is their Poisson bracket 
{/, 9}, as follows from the Jacobi identity, {{/, g}, 1} = {/, {g, I}} + {g, {I, /}} = 0. 
Thus, this Poisson bracket can be computed directly in angular momentum space 
without going back to the bracket (j23|) ; the result is the Lie-Poisson bracket, 



Interpretations of these spaces may be given in terms of the theory of "reduction" 
(Marsden and Ratiu, 1999). Angular momentum space is the Poisson manifold that 
results from Poisson reduction of the phase space R 4 under the U(l) action (|27|) 
generated by /. It is not by itself an ordinary phase space (symplectic manifold), which 
would have an even dimensionality, but it is foliated into symplectic submanifolds 
(the symplectic leaves) . In this case the symplectic leaves are the 2-spheres in angular 
momentum space, that is, the level sets J 2 = j 2 , the images under tt of the 3-spheres 
I = j in R 4 or C 2 . Canonical coordinates on a given 2-sphere J 2 = j 2 are (<fi,J z ), 
a (q,p) pair, where J z — jcosQ and where (9, </>) are the usual spherical angles in 
angular momentum space. Thus we have 



and the symplectic form on a given sphere is j c?f2, where dQ is the element of 
solid angle. This is not the geometrical solid angle in a Euclidean geometry on 
angular momentum space, which would be j 2 d£l. Another interpretation of angular 
momentum space is that it is the dual of the Lie algebra of SU(2), while tt, given by 
(|2"5j) . is the momentum map of the SU{2) action (f2"9"|) . 

We now have three spaces, the "large phase space" K 4 or C 2 , its image under 
7r, "angular momentum space" K 3 , and its symplectic leaves, the "small phase 
spaces," the 2-spheres J 2 = j 2 . Angular momentum space is useful for visualizing 
classical angular momentum vectors, but by considering inverse projections under tt 
the corresponding geometrical objects in the large phase space can be constructed. 




(30) 



dq A dp = d<f) A d(j cos 9) — j sin 9 d9 A dip — j d£l, 



(31) 



Semiclassical analysis of Wigner 3j -symbol 



12 



Angular momentum space has been used since the time of the old quantum theory for 
visualizing the classical limit of quantum angular momentum operators; for example, 
one spoke of an angular momentum vector "precessing" around the z-direction. In 
reality, the "precession" defines a manifold of classical states in the small phase space 
that is a level set of a complete set of commuting observables, that is, it is an invariant 
torus of an integrable system (just a circle in the lj-model, where the commuting 
observables are I and J z ). 



5.2. The Nj -model 

We now consider the classical mechanics of the iVj-model, which is mostly a simple 
generalization of the lj-model. We have 2N classical oscillators H ril = (l/2)(x 2 + 
Prfj); the configuration space is (R 2 )^ = M. 2N and the "large" phase space is (M 4 )^ = 
R 4N or (C 2 ) N = C 2N . We define z r ^ = (x r „ + ip r „)/V2, z r „ = (x r „ - i Pr ^/V2, 
so a point in phase space can be thought of as a collection of iV 2-spinors, (z r i, z r2 ), 
r = l,...,N. We make the obvious definitions (classical versions of equations fl5])- 

to, 

H 

Jri 2 ^ ] C T \iv %rv: (33) 
[if 

as well as J 2 = J^i 4) Ji = J2 r J ™ or J = J2 r J r, and J 2 = J2i J t- 

We denote a value of the functions I r by j r > 0; for positive values j r > 0, 
r — 1, . . . , N, the level set I r — j r (or J 2 = j 2 ) in the phase space (M 4 )^ is 
S 3 x . . . x S 3 = (S 3 ) N . The flow generated by I r for a specific value of r is just 
multiplication of the r-th spinor (z rl , z r2 ) by a phase factor exp(— iip r /2), as in l[27p: 
the other spinors are not affected. Thus the N commuting flows generated by all the 
J r 's constitute a U(1) N = T N action on the large phase space (T N is the iV- torus). 

Equation (I551> defines the projection map ir : (C 2 ) N — > (R 3 )^, the latter space 
being "angular momentum space" for the Nj -model, with one copy of IR 3 for each 
classical angular momentum vector J r . In view of its importance, we write out the 
components of this map explicitly: 

Jrx = -^{ZrlZ r 2 + Z r2 Z rl ) = Re(z rl Z r2 ), (34) 

1 

J ry = -^{ZrlZr2 - Z r2 Z r i) = lm(z rl Z r2 ), (35) 
Jrz =i(|z rl | 2 -|z r2 | 2 ). (36) 

Points of angular momentum space can be visualized as N classical angular momentum 
vectors, each living in its own angular momentum space, or iV such vectors all in the 
same 3-dimensional angular momentum space. The inverse image under 7r of a set of 
N nonvanishing classical angular momentum vectors is an A^-torus in the large phase 
space, generated by taking any point in the inverse image (a collection of N 2-spinors), 
and multiplying them by N independent, overall phase factors. We denote the angles 
on this torus by tp r , r — I, . . . , N, which are the evolution parameters corresponding 
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to the I r , as in (|27[k thus their periods are 47r. As in the lj-model, angular momentum 
space (R 3 )^ is a Poisson manifold, now with Poisson bracket 

The symplectic leaves (the "small phase spaces") are the spaces S 2 x . . . x S 2 — (S 2 ) N 
obtained by fixing the values of ji, ■ ■ ■ , j'jv, with canonical coordinates ((f> r , J rz ) on 
each sphere. 

In the classical iVj-model any partial or total sum of the angular momenta 3 r 
generates an SU(2) action on the large phase space, generalizing equations ([28]) and 
(129|) in the lj-model, in that the ST/ (2) matrix u is applied to all spinors (z r i,z r 2) 
whose r values lie in the sum. For example, the total J rotates all spinors. 

These SU(2) actions on the large phase space project to SO(3) actions on angular 
momentum space. Consider, for example, the SU(2) action generated by the total J. 
Along an orbit in the large phase space generated by n • J, parameterized by 9, we 
can follow the value of J r , giving us 3 r (ff), an orbit in the small phase space (the 
projection under 7r of the first orbit). The latter orbit is 

3 

where R(n,6) is the 3x3 rotation associated with u(n, 6) according to (f2"Tj) . This 
is the classical analog of (|20p . It follows from (|33[) and the spinor adjoint equation, 
iiV'u = ^2 j Rij o-i , itself equivalent to (f2"Tj) . Thus, under the SU(2) action on the 
large phase space generated by J, the individual vectors 3 r rotate in their individual 
angular momentum spaces by the corresponding 3x3 rotation. For example, J z 
rotates all vectors J r about the z-axis. Because of the two-to-one relation between 
SU(2) and SO(3), when the orbit in the large phase space goes around once (9 goes 
from to 4tt), the angular momentum vectors go around twice in their individual 
angular momentum spaces. 

We may visualize this action as in Fig. [TJ where A represents a point of angular 
momentum space (a set of N classical angular momentum vectors 3 r in the Nj- 
model). To obtain the generic case we assume these vectors are linearly independent 
(in particular, none of them vanishes). In the figure, T is the inverse image of A under 
7r, an TV-torus. Point a is any specific point in the large phase space on this iV-torus, 
to which the SU{2) rotation u(n, 9) is applied for < 9 < 4ir. That is, we treat a as 
initial conditions for the Hamiltonian flow generated by n-J, with 9 as the parameter. 
This generates the circle C in the large phase space, which amounts to rotating all N 
spinors by the same u(n,6). The projection of the circle C is a circle D in angular 
momentum space (IR 3 )^, that is, all classical vectors 3 r rotate about n by angle 9. 
However, when the circle C is covered once, circle D is covered twice. This is because 
when 9 = 2-7T, the spinor rotation u(n, 9) = —1, so all spinors in the large phase space 
are just multiplied by —1. This is illustrated as point a' in the figure, where all spinors 
are —1 times their values at a. Since —1 is just a phase factor, both a and a' project 
onto the same point A in angular momentum space. These are the only two points on 
C that project onto A; for 9 not a multiple of 2ir, the spinor rotation u(n, 9) is not a 
multiplication by a phase factor. 

Alternatively, we may apply the entire group SU{2) to the original point a (not 
just rotations along a fixed axis). Then the manifold C is the orbit of the SU(2) 
action which is diffeomorphic to SU{2). The point a 1 is the image of a under u = — 1, 
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a specific element in SU(2), and once again it projects onto the original point A 
in the small phase space. The manifold D is the orbit of point A under the group 
SO(3). In the lj-model, it is just a sphere in angular momentum space (all vectors 
that can be reached from the original one by applying all rotations) , while in the Nj- 
model for N > 1 D is generically diffeomorphic to SO(3) (it is the set of all classical 
configurations of N angular momentum vectors that can be obtained from the original 
one by applying rigid rotations). 

6. The invariant jm-tori 

In this section we continue with the classical point of view, examining the classical 
manifolds corresponding to the left side of the matrix element (f2"2")l . For this exercise 
and the rest of the paper we adopt a 3j-model (N = 3). The manifolds in question 
are the level sets of the commuting functions I r , J rz , r = 1,2,3, or, equally well, 
of the functions J rA , = (l/2)|z rA1 | 2 for r = 1,2,3, fj, = 1,2, since I r — I r i + I r 2 and 
Jrz — Iri — Ir2- We denote the level sets by I r = j r , J rz = m r for contour values j r , 
m r , r = 1,2,3, or, equivalently, by 

Iri = 7j(ir + m r), Ir2 = ~(jr - «V)' (39) 

In spite of the notation, j r and m r take on continuous values and are not necessarily 
the eigenvalues of any quantum operators. Since the I rfl are all nonnegative, we have 

jr > 0, -j r < m r < jr, (40) 

the classical analogs of the usual inequalities in quantum mechanics. 

Since each of the six I rji is a harmonic oscillator (times 1/2), the level set of 
the .Zrp's is an invariant torus of a collection of harmonic oscillators. Generically (for 
nonzero amplitude in each oscillator, that is, when none of the quantities j r ± m r 
vanishes) this is a 6-torus, upon which the coordinates may be taken to be the six 
angles 9 rfi , the variables of evolution of the I r ^. The Hamiltonian flow generated by 
I r{i for a specific value of r and /j, just multiplies z r ^ for the same values of r and /i 
by exp(— i6 r ^/2) 1 while leaving all other z's unaffected. This is not an overall spinor 
rotation since the other half of the spinor containing the given z rfl is not affected. 
If viewed in the Cartesian x rfl -p r ^ phase plane, this flow is a clockwise rotation by 
angle 9 rll /2, as illustrated in Fig. O The period of the angles 9 rfjL is An. We agree to 
measure the angles 9 rfl from the positive x ril axis, as in the figure, where z rfl is real 
and positive (or zero); this is a specific convention for a set of canonical coordinates 
(9 r ^,Ir^), T = 1, 2, 3, [i = 1, 2 on the large phase space. The volume of the 6-torus 
with respect to the measure d9u A ... A d9 32 is (47r) 6 . 

These tori are also the orbits of the flows generated by the observables I r , J rz . 
We denote the evolution variables of the I r and the J rz by tp r (as above) and <p r , 
respectively. Each J rz generates an SU{2) rotation about the z-axis on the spinor 
with the given value of r; thus each of the six angles (rp r , <p r ) has period An. However, 
when we allow all six angles (?p r ,(f) r ) to range from to An, the torus is actually 
covered eight times. This can be seen from Fig. [TJ a rotation by 2n in one of the (f>'s 
and one of the "0's returns us to the initial point (the path is a to a' along C in Fig.[TJ 
then a 1 to a along T.) Alternatively, we may consider the canonical transformation 
(9 rf *; I rl i) -> (Vv, (f> r ;I r , Jrz), generated by 

F 2 (9 rl ,9 r2 ,I r ,J rz ) = hfi r l{I T + Jrz) + 9 r2 (I r - J rz )}, (41) 
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which generates I r = I r \ + I r 2, J, 



rz 



I r \ — I r 2 and 

(fir = 7;(6rl — Qr2), 



(42) 



so that the Jacobian in the angles is (1/2) 3 = 1/8. To cover the torus precisely once 
we may let the Vv's range from to 47r and the r 's from to 2-7T, or vice versa; thus 
the volume of the torus with respect to difii A dip2 A dips A d(fi\ A d(j>2 A dfo is 



The angles <fi r defined in this way on the large phase space can be projected onto 
the small phase space, whereupon they coincide with the usual azimuthal spherical 
angles in the individual angular momentum spaces. It is clear this must be so to within 
an additive constant, since the J r2 -flow is just an S*0(3) rotation about the z-axis in 
the r-th angular momentum space, but by our conventions even the additive constant 
comes out right. To see this we note first of all that (fi r for a given r is constant 
along the Inflows for all s, since the variables in question are members of a canonical 
coordinate system on the large phase space and satisfy {(fi ri I s } = 0. Thus, (f> r , defined 
in the large phase space, projects onto a meaningful function in angular momentum 
space. Next, to compute the value of (fi r for a specific angular momentum vector J r , 
it suffices to take any point in the 3-torus that is the inverse image, that is, any value 
of the angles ip s may be chosen. For simplicity we take ifi r — 0, which implies 9 r i = <fi r 
and 9 r 2 — —(fir- This in turn implies z r \ = \z r \\ exp(— i(fi r /2), z r 2 = \z r 2\ exp(i0 r /2). 
But by equations ([34]) and (|35|) . these imply J rx — J r ± cos</v, J ry — J r j_ sin</v, where 

Jr± = \Zrl\\z r 2\- 

We shall henceforth call the level set I r = j r , J rz = m r the "jm-torus." This torus 
can be projected onto angular momentum space; we consider the generic case when 
j r ± m r for all r, in which case the jm-torus is a 6-torus. In this case, its image in 
angular momentum space is a 3-torus, which, since it is a surface on which I r — j r , is 
also a submanifold of the small phase space. This is because the three ip r angles just 
change the overall phases of the three spinors, without changing their image under n, 
so the three coordinates on the projected 3-torus are the angles (fir- The 3-torus in 
angular momentum space can be visualized as three classical vectors J r in a single 
angular momentum space, with specified values of m r = J rz , "precessing" about the 
z-axis. See Fig. [3] This is an example of how we shall visualize manifolds in the large 
phase space: The jm-torus, a six-dimensional manifold in the large phase space (itself 
with twelve dimensions), is visualized as three angular momentum vectors in three 
dimensional space, as in Fig. [3l defining a 3-torus by varying their azimuthal angles 
independently, and each point of this 3-torus is associated with another 3-torus, the 
inverse projection under n of the given point, which consists of independently changing 
the overall phases of the three spinors. The 6-dimensional jm-torus is thus conceived 
of as the Cartesian product T 3 x T 3 (it is actually a trivial bundle). 

7. The Wigner manifold 

Now we turn to the right hand side of the matrix element (|22|) . containing the state 
IJ1J2J3O). This state suggests that we examine the classical manifold upon which the 
I r have definite values, say, I r = j r , and upon which J 2 = 0. Again, we do not 
necessarily identify the j r with any quantum numbers, but it is convenient in the 
following to assume that none of the j r 's vanishes. 



V jm = (2tt) 3 {Att) 3 . 



(43) 
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7.1. Properties of the Wigner manifold 

Usually the dimensionality of a manifold can be guessed by counting the constraints 
that define it, for example, we expect the manifold in the large (twelve-dimensional) 
phase space upon which I r = j r , Jf 2 = j 2 2 , J z = m and J 2 = j 2 , for given contour 
values, to be six-dimensional (six constraints on twelve variables). Indeed, for most 
values of j this is correct, and the manifold in question is a 6-torus (by the Liouvillc- 
Arnold theorem, for certain ranges of the contour values). These are the invariant tori 
that would be involved in the semiclassical treatment of the addition of three angular 
momenta, producing a nonzero result (the case j ^ 0). But this naive dimension count 
only works when the differentials of the functions in question are linearly independent 
(in particular, nonvanishing) on the manifold. This condition breaks down when 
J 2 = j 2 = 0, since J 2 = and J = imply one another, and d(3 2 ) = 23 ■ dJ = 0. 
In fact, just the four conditions I r = j r > 0, J 2 = j 2 = define a six-dimensional 
manifold in the twelve-dimensional large phase space (for certain ranges of the j r ). To 
see this we notice first that since J 2 = implies Ji = 0, i = 1,2,3, and J z = m = in 
particular, the J z constraint is not independent; neither is the J 2 2 = j 2 2 constraint, 
since when J = 0, j\2 = J3- This is just as in the quantum case. In fact, the manifold 
Ir = jr, J 2 = is characterized cquivalcntly but better by I r = j r , J = 0, since the six 
differentials dl r and dJi are linearly independent on it. (Although the Jj vanish on the 
manifold in question, their differentials do not.) Thus, the naive count of dimensions 
works with the set I r , Ji. 

We shall call the manifold I r — j r , J = in the large phase space the 
"Wigner manifold," because it corresponds to the rotationally invariant state | ji J2J3O) 
introduced by Wigner in his definition of the 3j-symbols. The dimensionality of this 
manifold (six, in the appropriate ranges of the j r 's) is the same as that of the invariant 
tori of any integrable system of six degrees of freedom, and indeed the same as that of 
the nearby invariant tori in phase space corresponding to the level sets of the functions 
(Ii, I2, 13, 3 12, Jz, J 2 ) when j > 0. The Wigner manifold, however, is not a torus. This 
is not a contradiction of the Liouvillc- Arnold theorem, which requires that the classical 
observables making up the level set should Poisson commute. In the case of the Wigner 
manifold, we do have {I r ,I s } = and {I r , J{\ = 0, but {Ji,Jj} — ~^2 k ^ijk Jk- The 
Hamiltonian vector fields corresponding to the functions I r , Ji are linearly independent 
on the Wigner manifold, but the three J^-flows do not commute (two Hamiltonian flows 
commute if and only if their Poisson bracket is a constant). The Wigner manifold is, 
however, an orbit of the collective action of these Hamiltonian flows (any point can 
be reached from any other point by following the flows in some order). These facts, 
information about the topology of the Wigner manifold, and the required ranges on 
the contour values j r will be clarified momentarily. 

The Wigner manifold is also a Lagrangian manifold, like the invariant tori of an 
integrable system. This means that the integral of p dx along the manifold is locally 
independent of path, so an action function S(x, A) can be defined. This function in 
turn is the solution of the simultaneous Hamilton- Jacobi equations for the observables 
(Ji, I2, h, J x , Jy, Jz), call them A,, i = 1, . . . , 6 for short, of which the Wigner manifold 
is the level set. 

To prove that the Wigner manifold is Lagrangian, we note that the differentials 
dAi are linearly independent, so the vector fields Xi are, too, and span the tangent 
space to the Wigner manifold at each point. Evaluating the symplectic form on these 
vector fields, we have u)(Xi,Xj) = — {Ai,Aj}. These Poisson brackets all vanish 
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except for the {Ji,Jj}; the latter are nonzero at most points in phase space, but 
on the Wigner manifold where J = 0, these also vanish. Thus the symplectic form 
restricted to the Wigner manifold vanishes, the condition that the Wigner manifold 
be Lagrangian. 

To visualize the Wigner manifold we work our way up from angular momentum 
space to the large phase space. First we attempt to construct three angular momentum 
vectors of given positive lengths ji, j'3 that add up to the zero vector. This can be 
done if and only if the j r satisfy the triangle inequalities, whereupon the values of the 
jV's (the lengths of the sides) specify a triangle that is unique to within orientation. If 
we choose a standard or reference orientation for the triangle, then the three desired 
vectors are the vectors running along its sides. Let us therefore assume the triangle 
inequalities are satisfied, and let us choose a standard orientation for the triangle by 
placing the J3 along the z-axis, Ji in the x-z plane with J\ x > 0, and J2 in the x-z 
plane with J^ x < 0, as illustrated in Fig. HJ Given any two triangles with the same 
(positive) sides, there exists a unique rotation that maps one into the other; this fact 
and others regarding triangles are discussed in the context of the 3-body problem by 
Littlejohn and Reinsch (1995). In the present context this means that all classical 
configurations of three classical angular momentum vectors of fixed lengths that add 
up to the zero vector are related to any one such configuration, such as the one shown 
in Fig. SI by a unique rotation. Thus the manifold of such classical configurations in 
angular momentum space I 3 x I 3 x R 3 or in the small phase space S 2 x S 2 x S 2 is 
diffeomorphic to SO(3). 

The Wigner manifold in the large phase space is now the inverse projection under 
7r of this SO(3) manifold in angular momentum space. Since the inverse image of any 
given point of angular momentum space is a 3-torus in the large phase space (obtained 
by varying the overall phases of the three spinors), the Wigner manifold is a 3-torus 
bundle over SO(3), and is six-dimensional. The bundle is nontrivial. 

The Wigner manifold may also be visualized with the help of Fig.0 an elaboration 
of Fig.[TJ It is assumed that the three j's are positive and satisfy the triangle inequality. 
The lower part of this figure refers to angular momentum space, while the upper part 
refers to the large phase space. Projection n maps between the two spaces. Point A 
in angular momentum space is a state of three classical angular momenta of the given 
lengths j r whose vector sum is zero (that is, the angular momenta define a triangle), 
in a definite orientation. To be specific, let us say that A is the configuration shown 
in Fig. |U By applying all SO(3) rotations to A we generate all orientations of the 
triangle, of which B in the figure is one. The lower circle in the figure represents the 
manifold of such configurations, diffeomorphic to SO (3). 

The inverse image of any point on this manifold under tt is a 3-torus in the large 
phase space. The 3-tori above points A and B are indicated schematically as lines 
Ta and Tb in the figure. Let a be some point on Ta- To be specific, if A is the 
configuration of angular momentum vectors shown in Fig. U and the stated conditions 
on the j r hold, then none of the three vectors lies on the negative z-axis. This means 
that for any point on Ta, \z r \\ 2 is never zero, since by equations 1|32[) and ([36|) . we 
have 

\Zrl? = jr + Jrz, \z r 2? = jr - Jrz , (44) 

for r = 1, 2, 3. Thus by adjusting the overall phases of the three spinors (z T \, 2V2), we 
can make z r \ real and positive for all r = 1, 2, 3. Let this be the point a in Fig. [5] 
It is notationally tempting to write m r for the value of J rz , but we shall not do 
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this in the context of the Wigner manifold, instead reserving the symbol m r for the 
contour value of J rz on the jm-torus. 

Now we apply spinor rotations to point a, that is, simultaneous multiplication of 
all three spinors (z r i, z r 2) by the same element of SU(2). The orbit thereby generated 
is a manifold diffeomorphic to SU(2), as indicated in Fig. [5l The projection of this 
manifold onto angular momentum space is the surface 50(3) shown in the figure, 
that is, all orientations of the triangle are generated. For example, in the 3-torus Tg 
over the angular momentum triangle with orientation B, there is a point b that can 
be reached from the given point a by some spinor rotation. The spinor rotation in 
question is one of the two that projects onto the SO(3) rotation that maps A into B, 
according to (|21|) . The orbit of reference point a under the SU{2) action therefore 
passes through the 3-tori over every possible orientation of the triangle. In fact, it 
passes through each 3-torus twice, since the SU{2) rotation u = — 1 is just a phase 
factor. This is the meaning of points a' and b' in the figure, which are related to points 
a and b by multiplying all three spinors by — 1. 

Thus any point on the Wigner manifold can be reached from the reference point 
a by applying some SU (2) rotation, and then adjusting the overall phases of the three 
spinors (z r x, z r 2). The first step is equivalent to following along the Hamiltonian flows 
in the large phase space generated by the three Ji (this creates the SU{2) rotation), 
while the second is equivalent to following the Hamiltonian flows generated by the three 
I r (this changes the overall phases of the three spinors). By letting the rotations range 
over all of SU (2) and the three angles ip r range from to 4ir, the Wigner manifold is 
covered twice. Thus we obtain coordinates on the Wigner manifold (a, j3, 7, ipx , ip2, ^3) 
(the first three of which are Euler angles on SU(2)). 

Solving the simulataneous amplitude transport equations for the six observables 
I r , Ji, r,i = 1,2,3 requires us to find an invariant measure on the Wigner manifold, 
that is, one invariant under all of the corresponding Hamiltonian flows. Details are 
presented in Sec. [Til for now we just guess that this measure is the Haar measure on 
the group SU(2) times the obvious measure on the 3-tori generated by the I r , namely, 



where (a,/3, 7) are Euler angles on SU(2). The integral of this measure over the 
Wigner manifold is 



where the 1/2 compensates for the fact that the Wigner manifold is covered twice 
when the Euler angles run over SU(2) and each ip r runs from to 4-7T. 

7.2. Angles related to the shape of the triangle 

Figure H] defines the angles rj r as the angles opposite vectors J r . Under our 
assumptions, these angles lie in the range < r\ r < it. By projecting all three vectors 
onto the directions parallel and orthogonal to each of the vectors in turn, we obtain a 
series of identities, 



sin 13 da A df3 A A dtpi A dip 2 A dip3 , 



(45) 



TV = -(167r 2 )(4^) 3 = 2V, 



(46) 



ji cos 772 + n cos 771 + j 3 = 0, 
ji sin 772 - j 2 sin = 0, 



(47) 
(48) 
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and four more obtained by cycling indices 1,2,3. These allow us to solve for the 
cosines of the angles rj r , 

C0S7?1 = -I" ->\ (49) 

2j2j3 

and cyclic permutations, which in view of the stated ranges on the angles allows all 
three angles r\ r to be uniquely determined as functions of (ji, 321 3%) ■ We shall regard 
angles rj r as convenient substitutions for these definite functions of the lengths of the 
angular momentum vectors. 

For the stated ranges on the r] r , the sines of the angles are nonnegative and are 
related to the area A of the triangle, as follows: 

A = — j Ji x J 2 | = ^jij2smn 3 

= jVih + h + .h) + 32 + 33) (ji - 32 + + 32 - h), (50) 

and cyclic permutations. Some authors define A as the final square root (without the 
1/4). 



8. Intersections of manifolds 



The stationary phase points of the matrix element (|22|) are the intersections of the 
jm-manifold and the Wigner manifold in the large phase space. Thus we must use a 
version of ([5]) for the matrix element, rather than ([3]). In this section we study the 
intersections of the manifolds, continuing with a classical picture. 

If the jm-torus and the Wigner manifold in the large phase space have a 
common point of intersection, then the projections of these two manifolds onto angular 
momentum space must have a common point of intersection. The converse is also true: 
if the projections have a point in common, then the inverse image of this point under 
7r, a 3-torus which is the orbit of the three 7,,-flows, must contain two points, one of 
which belongs to the jm-torus, and the other to the Wigner manifold. But the 3-torus 
is the orbit of the / r -fiows, and these flows confine one to both the jm-torus and the 
Wigner manifold. Therefore the entire 3-torus is common to both the jm-torus and 
the Wigner manifold. Therefore to find intersections of the jm-torus and the Wigner 
manifold, we may first find the intersections of their projections under tt. 

The jm-torus projects onto a set of configurations of three angular momenta 
with given lengths and fixed values of m r = J rzi with arbitrary azimuthal angles, 
while the Wigner manifold projects onto configurations with the same lengths in which 
the vector sum of the angular momenta vanishes, forming a triangle, with arbitrary 
orientation. Therefore to find a common point between these two sets of classical 
angular momenta configurations, we can either adjust the azimuthal angles of the 
angular momenta with given m values until a triangle is formed (total J = 0), or we 
can rotate a triangle from a given, reference orientation until the m values are the 
desired ones. We choose the latter procedure. 

Our reference orientation of the triangle is shown in Fig. [4l which is indicated 
schematically as the point A in Fig. [5j We must rotate this reference orientation to 
obtain some prescribed values of m r . These values satisfy the relations (|40[) . so in 
particular |m,3| < j 3 . Thus by rotating the triangle in the reference orientation about 
the y-axis by a unique angle f3, < (3 < tt, defined by 

m 3 =j 3 cos/3, (51) 
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we guarantee that J 3 has the right projection. The result of this rotation is shown in 
Fig. [51 for a certain negative value of m 3 . 

Next we rotate the triangle about the axis J3 by an angle 7, which does not 
change J3 or its projection, but which rotates Ji and J2 in a cone, as illustrated in 
Fig. [7j We wish to choose the angle 7 so that J 2 has the desired projection 777,2 onto 
the 2-axis. In Fig. [7j Ji is not shown, but J2 rotates about the J3 direction, its tip 
sweeping out circle C3. The circle C z in the figure is swept out by a vector of length 
32 and projection 777,2 onto the z-axis (this vector is not shown). Circles C3 and C z 
intersect in two points Q and Q' in the figure, which represent two orientations of the 
triangle that have the correct values of both 7713 and 777,2 . Now the orientation of the 
triangle is fixed, so there is no more freedom to rotate J\, In this final orientation 
the value of J\ z is — J22 — Jzz = —m2 — ma, since J = for the triangle. Either 
this value of J\ z equals the value of mi associated with the jm-torus, or it does not. 
If it does not, then there are no intersections between the jm-torus and the Wigner 
manifold. This is just the classical expression of the condition that the matrix element 
(|22|l vanishes unless m r — 0- Henceforth we assume that the m r values for the 
jm-torus do satisfy this condition. 

In this case we may solve for the values of 7 associated with points Q and Q' in 
Fig. [7J Writing i?(n, 9) for a 3-dimensional rotation by angle 9 about axis n, we have 
applied the rotation 

flOia. 7)fl(y, P) = R(y, P)R(*, i) (52) 

to the reference orientation in Fig. 21 where J3 is the unit vector in the direction J3 
shown in Fig. [5] (after the first rotation) . In the reference orientation the vectors are 

(sinf?2 \ / -sin 771 \ / 

, J 2 =ja , J 3 =j 3 I I . (,>;-!) 
cosry 2 / \ cos?7i / 

After applying rotation (|52)l these become 

cos P cos 7 sin 772 + sin /3 cos r/2 
Ji = ji I sin 7 sin 772 

— sin (3 cos 7 sin 772 + cos (3 cos 772 

— cos j3 cos 7 sin 771 + sin f3 cos 771 
1 2 = ,1 2 \ - sin 7 sin 771 

sin j3 cos 7 sin 771 + cos (3 cos 771 

sin/3 \ 

! . (54) 
cos (3 J 

We have already solved for (3 in (|5ip : we may now solve for 7 by demanding either 
J\ z = mi or J2 Z = 777,2- These lead to 

ji cos (3 cos 772 — mi 777,2 — 32 cos (3 cos 771 . 

cos 7 = , . . = . . . . 55) 

31 sin [3 sin 772 32 sin p sm 771 

These two conditions are equivalent (under the assumption m r — 0) > as follows 
from the identities (|47j) - (|49p . If the common value of the two expressions on the right 
hand side of (|55| lies in the range (—1, +1), then there are two real angles 7 satisfying 
(|55p . corresponding to the two points Q and Q' in Fig.[7j In this case the two manifolds 
have real intersections, and we are in the classically allowed region for the 3j-symbol. 
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We let 7 represent the root (the "principal branch") in the range [0, 7r], and —7 the 
root (the "secondary branch") in the range [— 71", 0]. Note that sin 7 > (< 0) on 
the principal (secondary) branch. If the right hand side of ((55)) lies outside the range 
[—1,1], then there are two complex roots for 7. In this case the two manifolds have no 
real intersections, but they do have complex ones. Only one of the two complex roots 
is picked up by the contour of integration used in obtaining the matrix element (|3]) or 
(|5j). resulting in an exponentially decaying expression for the matrix element. In this 
case we are in the classically forbidden region of the 3j-symbol. In the following for 
simplicity we assume we are in the classically allowed region. 

The points Q and Q' in Fig. [7] represent values of J 2 in a single angular momentum 
space R 3 . Taken with the values of Ji and J3, they specify points, call them P and P' , 
in the combined angular momentum space I 3 x I 3 x R 3 that lie on the intersection 
of the projections of the jm-torus and the Wigner manifold onto that space. Then 
by applying rotations about the z-axis to P and P', we generate a pair of circles in 
angular momentum space. Such rotations change neither the z-projections of the three 
vectors nor their vector sum (zero) . This is obviously a reflection of the fact that the 
operator J z defining the state on the right side of (f22|) is a function of the operators 
(Jiz, J2z, Jnz) defining the state on the left. Thus the projections of the jm-manifold 
and Wigner manifold under 7r intersect generically in a pair of circles. 

Thus the intersection of the jTO-torus and the Wigner manifold in the large phase 
space is the inverse image of this pair of circles under 7r, generically a pair of 3-torus 
bundles over a circle. Since the 7 r -flows and the J z -flow commute, these bundles arc 
trivial, in fact each is a 4-torus, on which coordinates are {tp\, ip2> ^3> <fi), where (ft is 
the angle of evolution along the J z -flow. The volume of either one of the 4-tori with 
respect to the measure dtfti A dift 2 A dtps A deft is 

V i = \{Ak)\ (56) 

the factor of 1/2 being explained by Fig. [TJ 

The jm-torus and the Wigner manifold, both six-dimensional, intersect 
in a 4-torus because the the lists of functions defining the two manifolds, 
{I\,l2,h,Jiz,J2z,J , iz) and {Ii,l2,I i i,J x ,Jy,Jz), have three functions in common 
while J z in the second list is a function of (Jiz, J22, Jsz) in the first list. Below 
we will transform the functions to make both lists have explicitly four variables in 
common (see Eq. (|97l)). 

9. Action integrals 

Action integrals on the jm-torus and the Wigner manifold are needed for the phases 
in expresssions such as (J5j) - We only need the action function at some point on 
the intersection between the two manifolds, which gives us a lot of choice since that 
intersection is a 4-torus. We continue with a classical picture in this section. 

9.1. Choosing reference points 

Action integrals are defined relative to some initial or reference point on each manifold. 
For the jm-torus, a convenient point is the one where 9 r ^ = 0, r = 1, 2, 3, \l = 1, 2, 
that is, the point where each z ril is real and nonnegative, as explained in Sec. O 
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According to (|39p. the spinors at this reference point are given explicitly by 

(57) 



Zrl \ _ ( Vj 



Z r2 J \ y/j r — m r 

The projection of this point onto angular momentum space is a set of vectors J r , 
r = 1,2,3 of given lengths j r that lie in the x-z plane, with J rz = m r and J x > 0, as 
shown by (f34|) - (f36|) . Such vectors are illustrated in Fig.O 

As for the Wigner manifold, it is convenient to take the reference point to be 
point a in Fig. which is discussed in Sec. [Jj This point projects onto the standard 
orientation of the triangle, point A in Fig. [5j where the angular momentum vectors 
have the values shown in (|53|) . At the point a, z r \ is real and positive for all r, as 
explained in Sec. [71 For example, for r = 1 this assumption combined with (j44|) implies 
zn = y/ji + Jiz, which by (|53)) becomes zu = + cos 772) = y/^Ji cos 772/2. Then 

(|35p and J\ y — imply that Z12 is purely real, and (|34| allows us to solve for zyi in 
terms of J± x , given by (f53|) . producing finally z%2 — y/^ji sin 772. Proceeding similarly 
with the other two spinors r = 2, 3, we obtain the three spinors at the reference point 
a on the Wigner manifold, 

z n \ _ nrr- ( cos 772/2 \ / z 2 i \ _ /y— / COS?7i/2 

Z12 J ~ V Jl ^ sin 772/2 J ' ^ 22 J - ^ -sinr7i/2 

' ^ Wsd). (58) 



^32 J V V . 

Now to obtain a point common to both the jro-torus and the Wigner manifold, 
we apply the spinor rotation 

f e-^/ 2 cos /3/2 -e^/ 2 smp/2\ 
M(y ' Wz ' 7)= ( e-^ 2 sin/3/2 e*r/» cob/J/2 J (59) 

to the reference spinors (z r i, ^2), r = 1, 2, 3, in (|58p . where Euler angles /3 and 7 are 
defined by (|5"Tj) and (|55|) . We obtain either the principal branch or the secondary one 
by taking 7 > or 7 < 0, respectively. The spinor rotation (f5"9"]) induces the 3x3 
rotation on the angular momentum vectors shown in (|52[) . Thus we obtain the spinors 
at the common point of intersection between the jm-torus and the Wigner manifold, 

' e- 47 / 2 cos /3/2 cos 772/2 -e 17 / 2 sin /3/2 sin 772/2 ' 
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(60) 
(61) 

= ' v-'M i. (62) 



212 J \ e-^l" 1 sin (3/2 cos 772/2 + e^/ 2 cos (3/2 sin 772/2 

'-^l 2 cos /3/2 cos 771/2 + e^l 2 sin /3/2 sin 771/2 



222/ \ e~ J7/2 sin (3/2 cos 771/2 - e 17 / 2 cos (3/2 sin 771/2 

z 3i ^ - P -ii/K/9i:{ cos/? / 2 



z 32 ; V ^ V sin /3/ 2 

One can easily check using (|3"4")) - (|3l))) that these spinors project onto the angular 
momentum vectors in (1541) . 

9.2. Computing the actions 

In computing action integrals we use the identity, 

^ \ Pr^J. dx r ^ — 2 ^ ^ (^r/^ dz r ^ Z r ^dz r ^) -\- — G? •ErnPrji • (63) 
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The integral of the left hand side is the usual action one would need for wave functions 
ip{xn, . . . , £32), but it can be replaced by the integral of the first differential form on 
the right, for the following reason. First, the integral of the exact differential on the 
right contributes the difference in the function (1/2) J^ru x rnPrn between the initial 
and final points. But the final point is the common point of intersection between 
the jm-torus and the Wigner manifold, so this contribution cancels when we subtract 
actions as in ([5]). As for the initial points on the two manifolds, these have been chosen 
(see Eqs. ([57)1 and ([55)) ) so that all z rfl are purely real, or p r/1 = 0. Thus the function 
in question vanishes at the initial points. As for the integral of the first term on the 
right of (I63)) . it can be written 

S = Im J z rfJ , dz rfl . (64) 

rfi 

For the action on the jm-torus between initial point (|57| and final point (I60j) - 
(|62p. we follow a path consisting of flows of the functions I TlL = (l/2)|z rA1 | 2 taken one 
at a time by angles 6 rfJi . Along the J r(1 -flow we have dz rfl /d6 r ^ — (i/2)z rfl , so the 
contribution to S is 

Im J -\z rfl \ 2 d9 r ^ = I r ^6 r ^, (65) 

since I rfJl is constant along its own flow and since 6 r{i = at the reference point. Thus 
the total action between initial and final points on the jm-torus is 

Sjm = ^ Irfi6rfi- (66) 
rfi 

Under the canonical transformation (T4"Tj) this becomes 

S 3f = ^Zi 1 ^ + Jrz<Pr) = ^(jrA + m r (j) r ), (67) 
r r 

where in the final form we replace I r and J rz by their values on a given jm-torus. 

The angles 9 r ^ or (tp r ,(f> r ) are the coordinates of the final point specified by 
Eqs. (|5D)) - ([d^1) . The solutions of Hamilton's equations for the J r -flow can be written 
Zru,(@ru,) — z rfj,{0) exp(— i8 rfl /2) (see Fig. [5]) where the initial conditions are real and 
nonnegative, so we have 9 rfi = 2argz rAI . Combining this and (|42p . we can write the 
action on the jm-torus as 

Sj m = 2 ^2 Irn arg z rfl = j r aig(z rl z r2 ) + ^ m r arg(z r iz r2 ). (68) 

rfi r r 

Using Eqs. ([60|) -([62 |) . this can be written, 

Sjm = J37 + ji arg(cos (3 sin 772 + sin (3 cos 7 cos r]2 +i sin (3 sin 7) 

+ j2 ar g( — cos j3 sin rji + sin j3 cos 7 cos r\\ + i sin (3 sin 7) 

+ mi arg (sin (3 cos 772 + cos f3 cos 7 sin 772 + i sin 7 sin 772) 

+ 7712 arg(sin/3cos77i — cos j3 cos 7 sin 771 — i sin 7 sin 771). (69) 

Here we have used the rule arg(afr) = arg a + arg b, which is only valid for certain 
choices of branch of the arg function. A more careful analysis shows that is the 
correct action along a certain path from the initial to final point (the principal branch) 
on the jm-torus if the range of the arg function is taken to be [— tt, tt). (The path is 
defined by 7 < 611,622 < 27r, —7 < 612,621 < 7, that is, one integrates from to these 
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final 6 values.) In particular, this means that ipi , ip 2 7 "03 j 4>i au ne m [0, 7r], while 4>2 
lies in [— tt, 0]. These ranges on angles 4>i, 4>2 are also evident from Fig. [71 Similarly, 
for the secondary branch (—tt < 7 < 0, sin 7 < 0) there exists a path such that with 
the same range on the arg function ([M]) is still correct. With these understandings, 
the values of Sj m on the two branches differ by a sign. We shall henceforth write Sj m 
(~Sjm) for the principal (secondary) branch. 

Equation can also be written in terms of cos -1 functions. We note that 
arg(znzi 2 ) = cos _1 [Re(znZi 2 )/|2:iiZi2|] and that Iznz^l = (ji—nil) 1 ^ 2 , etc. We can 
also use (|55p to eliminate cos 7. For the principal branch (7 > 0) this gives 

_ . 1 / ?'i cos j3 777,1 cos 772 \ . _i / mi cos 771 — j 2 cos 8 

Sjm = Jl COS _ j ) + ]2 COS 1 



sinr7 2 Ji± / V sin77iJ 2 i 

s (3 cos t? 2 — mi \ _i I ' jl cos rj2 — it, 

■ ■ a ■ + mi cos ■ a t 

ji snip sin 7/2 / \ smpJi 

s r/i — rr 
sin /3J 2 



-1 i jl COS (i COS 7? 2 — 777,1 \ _i / jl COS 77 2 — 777,1 cos /3 

J3 COS ; — — + 7771 COS 

ji sin p sin 7/2 / 

1 (h COS 77i - 777 2 COS /3\ 

m 2 cos — — , (70) 



where 



J r ± = VJ? - (71) 

and where the range of the cos -1 function is [0, %]. Finally, by using Eqs. (|4D| and 
(|50p these can be written explicitly in terms of the parameters j ri m r . The result has 
the form of (|67|) . where 

= cos- 1 ( , (72) 
\ 4AJi^ / 

and cyclic permutations of indices, and where 

\ - jj ~ Jl ~ 277717773 \ ^ = _ cog _i / jl -j|-J 2 2 -2m 2 777 3 \ 
2Ji^J3J_ / ' V 2J2J_J3± / ' 

and where 03 = 0. 

Now we consider the action on the Wigner manifold between the initial point 
(|58p and the final point ([60)) (|62| . The path between these points is made up 
of the product of rotations ([55]) . so we consider the action integral (j64p along a 
rotation by angle 9 generated by n • J. Hamilton's equations (see Eq. (|28p) are 
dz ril /d6 = (i/2) J2„ z ru (n ■ a) Ufi , so by ffity we have 

r 9 i r e 

S = Im - y]z rv {n- a) Vfl z rfl d9 = / (n • J) d6 = (n • 3)9 = 0, (74) 

where we use (|33"1) . the fact that n • J is constant along its own flow, and the fact that 
J = on the Wigner manifold. The rotational action vanishes. 

Thus the phase of the matrix element (|22|) is determined entirely by the action 
integral along the jm-torus, that is, to within a sign it is given by Eqs. (1731) . 
This is the phase function determined previously by Ponzano and Regge, Miller, and 
others, and we see that it is essentially a simple combination of the phases of the 
Schwinger oscillators. We have, however, determined this phase function entirely 
within a classical model, that is, without imposing any quantization conditions on the 
manifolds. 
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10. Bohr-Sommerfeld quantization 

We do not need Bohr-Sommerfeld approximations to the eigenvalues of the operators 
involved in the 3j-symbols because those eigenvalues are known exactly. We must, 
however, quantize the jm-torus and the Wigner manifold, to obtain the wave functions 
whose scalar product is the 3j-symbol. We also need the Bohr-Sommerfeld rules to 
make the connection between the contour values for various classical functions and 
the standard quantum numbers of the associated operators. 

To quantize a Lagrangian manifold we must first find the generators of the 
fundamental or first homotopy group of the manifold, that is, a set of closed contours 
in terms of which all closed contours can be generated by concatening curves. In 
the following we shall call these generators "basis contours," although technically 
the fundamental group, even when Abelian, is a group and not a vector space. For 
example, in the familiar case of the invariant n-tori of integrable systems of n degrees of 
freedom, the fundamental group is Z", that is, an arbitrary closed contour is expressed 
as a "linear combination" of the n basis contours with integer coefficients. The n basis 
contours themselves go around the torus once in the n different directions. 

After finding the basis contours, we compute the total phase associated with 
each of them, the sum of an action, the integral of pdq around the contour, and a 
Maslov phase, which is —ir/2 times the Maslov index of the loop. Both these phases 
are topological invariants and are additive when loops are concatenated. Then we 
demand that the total phase be a multiple of 2ir; this is the consistency condition on 
the semiclassical wave function that selects out certain manifolds as being "quantized." 

The Lagrangian manifolds we are interested in are level sets of a set of classical 
functions that are the principal symbols of a set of operators, which in our application 
need not commute. Quantized Lagrangian manifolds support wave functions that are 
approximate cigcnfunctions of the set of operators. The corresponding eigenvalues are 
the contour values of the principal symbols, to within errors of order ft 2 . 

10.1. Quantizing the jm-tori 

In the case of the jm-tori, whose fundamental group is Z 6 , we are dealing with the 
eigenfunctions of a set of independent harmonic oscillators, so the problem could not 
be more elementary from a semiclassical standpoint. There are, however, interesting 
issues that arise. Let the complete set of commuting quantum observables be (I r , J rz ), 
r = 1,2,3, defined in ([5]) and Let us denote the Weyl symbol of an operator A by 
sym(^4). Then we have 

sym(/ r ) = I r - i sym(j r z) = J rz , (75) 

where I r and J rz (the classical functions) are defined by Eqs. (|32|) and (|33|) . The 
operator I r violates our assumption in Sec. [2] that the commuting operators defining 
our integrable system should have Weyl symbols that are even power series in ft 
(since the —1/2 is of order ft). Our assumption is valid for the harmonic oscillators 
H r = (1/2) Ylni^r/j, +Pr^i)> Du t m defining I r in (JSJ we have subtracted the zero point 
energy, a constant of order ft, I r = (1/2) (if r — 1), so that the eigenvalues of I r would 
be the conventional quantum numbers j r for an angular momentum, and so that the 
identity (fT3|) would have a familiar form. In the following we shall take the principal 
symbol of I r to be the whole symbol, including the —1/2. This achieves the same 
results we would have had if we had worked with H r instead of I r and defined the 
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principal symbol as the leading term in h, as in Sec. [2J since sym(_ff r ) = 2I r . We must 
be careful, however, since the principal symbol of I r is not I r . For most of the other 
operators we shall use, the principal symbol is obtained simply by removing the hat 
(for example, J rz above). 

The basis contours on the jm-torus are most easily expressed as the contours on 
which each of the angles 9 rf i is allowed to go from to Aw while all other 9 rfJ _'s are 
held fixed. We may also use any linear combination of these contours with integer 
coefficients and unit determinant. In terms of the angles Vv, 4>r, given by (|42|) . a 
convenient choice is to take one basis contour as the path on which one tp r goes from 
to 4-7T while all others ?/v's and all <^> r 's are held fixed; this is following the / r -flow for 
elapsed angle ijj r — Air. A second basis contour may be taken to be the path on which 
ip r goes from to 2tt, and then tj> r goes from to 27r; the two legs involve following 
the J r -flow and then the J rz -How, each for elapsed angle 2tt. Doing this for r = 1, 2, 3 
gives us six basis contours on the jm-torus. 

The action along the first basis contour is computed as in (f6"5)) . Hamilton's 
equations for I r are dz rtl /dip r = (i/2)z rM , so we obtain S r i = 47rJ n where S r \ refers 
to the action along the first basis contour. For the second basis contour, the first leg 
contributes an action 27r/ r , while the second leg, which follows the flow generated by 
J rz , is a rotation whose action may be computed as in (|74|) . but with J replaced by J r 
since we do not sum over r. Thus the final answer does not vanish (J r is nonzero on the 
jm-torus), and the contribution from the second leg is 2i:J rz , or S r 2 = 2n(I r + J rz ). 
Altogether, we have 

S r i = 47rj r , S r2 = 2n(j r + m r ), (76) 

where the 1 and 2 refer to the first and second basis contours associated with a 
particular value of r, and where we have replaced I r and J rz by their contour values 
j r and m r on the jm-torus. 

Next we need the Maslov indices along the two basis contours. Here we follow 
the computational method described in Littlejohn and Robbins (1987), which uses 
the determinant of complex matrices and which is based ultimately on Arnold (1967). 
Similar techniques are discussed by Mishchenko et al (1990). The method works for 
finding Maslov indices along closed curves on orientable Lagrangian manifolds in R 2 ™ . 
To describe the method we adopt a general notation, in which global coordinates on 
phase space are (q±, . . . , q n ,Pi, ■ ■ ■ ,p n )- We suppose that there exists a set of n vector 
fields on the Lagrangian manifold, linearly independent at each point, so that they 
span the Lagrangian tangent plane at each point. In our applications, these are the 
Hamiltonian vector fields associated with a set of functions (A\, . . . , A n ). We consider 
the rate of change of the quantities qi — ipi along the j-th vector field, which is the 
Poisson bracket {qi — ipi,Aj}. The set of these Poisson brackets forms an n x n 
complex matrix Mij that is never singular, so det M traces out a closed loop in the 
complex plane without passing through the origin when we go around a closed loop 
on the Lagrangian manifold. Then the Maslov index /i associated with this loop is 
given by 

/i = 2wndetM, (77) 

where wn refers to the winding number of the loop in the complex plane, reckoned as 
positive in the counterclockwise direction. The winding number is invariant when My 
is multiplied by any nonzero complex constant (or constant matrix), so such constants 
can be dropped in the calculation. 
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For the jm-torus, we identify the g's and p's with the coordinates x rfl and p rfJi , 
and the A's with the functions (7i, I21 13, Jiz-, J2z-, Jzz)- Then we can replace ^ — ipi 
by z rfJ/ , dropping the 1/V2. The needed matrix elements are 

OZ r ^ Z 

{Zrl, J sz } = = -5 rs Z r l, ( 78 ) 

a z r \ 2 

S- 7 1 _ - dJsz - i X ~ 
OZ r 2 2 

We drop the constant i/2 on the right hand side, and choose the ordering 
(h, Jiz, I2, J22, I3, J3z) for the functions. Then the 6x6 matrix block diagonalizes 
into three 2x2 blocks, and we find 

detM = Z11Z12Z21Z22Z31Z32, (79) 

to within a constant. 

Along the flow of one of the I r 's we have Z r u,(ipr) = ex p(*'0r-/2)zr/j(O), ° r 
detM(^v) = exp(ii/v) det M(0). Therefore when the given ip r goes from to 47r, 
the other f/v's being held fixed (this is the first basis contour), detM circles the 
origin twice and we have /i = 4. Along the flow of one of the J rz 's, however, we 
have z r i(<pr) — exp(i0 r /2),? r i(O) and z r 2{4> r ) = exp(— i0 r /2)z r 2(O), or detM(<ft r ) = 
detM(O). Therefore along the second basis contour the J r -flow takes detM once 
around the origin (elapsed parameter ip r = 27r), while the J rz -&ow does nothing. 
Therefore the Maslov index of the second basis contour is /i = 2. There are easier 
ways to find the Maslov indices of harmonic oscillators, but this calculation is useful 
practice for the case of the Wigner manifold that we take up momentarily. 

Now we apply the quantization conditions. For the first basis contour the total 
phase is 4irj r — 4(7r/2), which we set to 2n r n where n r is an integer. Thus the quantized 
tori must satisfy j r — (n r + l)/2. The allowed values of n r are determined by the fact 
n r < — 1 is impossible in view of the fact that I r is nonnegative definite, and n r = — 1 
corresponds to a torus of less than full dimensionality (six) , so the wave function (fTJ) 
is not meaningful. Thus we must have n r — 0, 1, . . .. The j r in these formulas, and 
throughout all of the classical analysis from Sec. [5] up to this point, has referred to a 
contour value for the function J r ; the only difference now is that we are restricting the 
value of j r in order that the torus be quantized. This j r , however, is not value of the 
principal symbol of the operator I r (see Eq. (1751 ). so the Bohr-Sommerfeld or EBK 
quantization rule gives the semiclassical eigenvalue of I r , call it as 

3T=jr-\ = \. (80) 

The semiclassical eigenvalues of I r are nonnegative integers or half-integers, the exact 
answer (not surprising in view of the fact that semiclassical quantization of quadratic 
Hamiltonians is exact) . If we use the operator identity (fT3)) to find the eigenvalues of 
operators J^, these are also exact. 

Equation (|80[) shows that the classical level set corresponding to quantum number 
is jr = j? 11 + 1/2- The extra 1/2 in this formula has caused some discussion in 
the past and merits a little more now. Ponzano and Regge (1968) used intuition 
and numerical evidence to argue for the presence of the 1/2. Miller, without 
knowing about Ponzano and Regge, also included the 1/2, referring to the "usual" 
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semiclassical replacement for angular momenta. Presumably he was referring to the 
similar replacement that occurs in the treatment of radial wave equations (the Langer 
modification, see Berry and Mount (1972), Morehead (1995)). It is not obvious to 
us what the Langer modification has to do with the 1/2 that occurs in the present 
context, nor are we aware of any general rules about when in the asymptotics of angular 
momentum theory it is correct to replace a classical j by j + 1/2 (instead of ' + 1)] 1 / 2 
or something else). Schulten and Gordon (1975b) and Reinsch and Morehead (1999) 
obtain the 1/2 as a part of their proper semiclassical analyses. Biedenharn and Louck 
(1981b) also speculate on the significance of the 1/2. Roberts' (1999) derivation of 
the asymptotics of the 6j-symbols does not produce the 1/2. He argues that in the 
asymptotic limit there is no confusion about whether a given point lies in the classically 
allowed or forbidden region, whether the 1/2 is included or not. The omission of the 
1/2 does, however, cause an error in the phase function that is of order unity, so 
the oscillations are not even approximately represented (this point was also made by 
Biedenharn and Louck, 1981b). We suspect that a modification of Roberts' method 
would produce the 1/2's. According to Girelli and Livine (2005), different choices for 
the semiclassical replacement for the quantum number j have been made by various 
researchers in the field of quantum gravity. Here we have shown that the extra 1/2 
is a necessary consequence of standard semiclassical theory. We remark in addition 
that with the inclusion of the 1/2, the quantized spheres in angular momentum space 
are those with an area of (2j qu + l)27r, that is, they contain a number of Planck cells 
exactly equal to the dimension of the irrep, obviously a form of geometric quantization. 
In particular, the s-wave j qu = is represented by a sphere of nonzero radius, a case 
for which the replacement j r — j qu + l/2 is declared by Biedenharn and Louck (1981b) 
to be "clearly invalid." 

For the second basis contour on the jm-torus, the quantization condition is 
27r(jV + rn r ) — 2(7r/2) = 2irn' r , where n' r is an integer. With (f80|) . this implies 
m r = — j' qu + n' r . Combined with the classical restriction (j40| . this gives the usual 
range on magnetic quantum numbers. Again the semiclassical quantization is exact. 
In the case of J rz , the eigenvalue of the operator is equal to the classical contour value 
m r on the quantized torus (without any correction such as we see in (f8"0"|) ). 

10.2. Quantizing the Wigner manifold 

We begin the quantization of the Wigner manifold by guessing the basis contours of 
the fundamental group by inspection of Fig[5j Taking the base (initial) point of the 
loops to be point a in the figure, we get three independent basis contours (call them 
Ci, C2, C3) by going around the 3-torus Ta in the three different directions. A fourth 
contour (call it C4) is created by following an SU(2) rotation about some axis by 
angle 2tt, taking us along the path aba' , which puts us half way around the torus Ta 
from the starting point, and then by applying half rotations along each of the three 
directions on the torus, taking us down along Ta in the diagram back to the starting 
point a. These four contours are not independent, since 

2d = Ci + C 2 + C 3 , (81) 
where addition of contours means concatenation, but they are convenient for studying 
the quanitzation conditions since a minimal set of three contours (not ((71,62,63), 
but for example (Ci, C2, C4)) is less symmetrical. The fundamental group is Z 3 . 

We may show the correctness of this guess by a topological argument. First, 
the Wigner manifold is the orbit of a SU{2) x T 3 group action on the large phase 
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space. Let {if>i,ip2,tp3) be coordinates on T 3 , where < ip r < Air. The action of 
(u, (ipi,ip2, V^)) G SU(2) x T 3 on the large phase space is to multiply each spinor by u 
and then by exp(— iip r /2). The isotropy subgroup of this action consists of the identity 
(1, (0, 0, 0)) and the element (—1, (2w, 2ir, 2%)), a normal subgroup isomorphic to Z2. 
Thus the Wigner manifold is diffeomorphic to (SU(2) x T 3 )/Z 2 , itself a group manifold, 
of which the original group SU{2) x T 3 is a double cover. This cover is topologically 
simple since 577(2) is simply connected (we can go to the universal cover if we wish 
by replacing T 3 by K 3 ). Therefore the homotopy classes on the Wigner manifold are 
in one-to-one correspondence with classes of topologically inequivalent curves that go 
from the identity in SU(2) x T 3 to one of the elements of the isotropy subgroup. Since 
SU(2) is simply connected, such paths are characterized by choice of the end point 
(the element of the isotropy subgroup), and the winding numbers around the torus T 3 . 
They are thus all "linear combinations" of the four contours Ck, k — 1, . . . , 4 defined 
above. Because of the relation (|5Tj) , however, the fundamental group is not Z 4 , but 
only Z 3 . 

It is easy to compute the action along these contours. Contours C r , r = 1,2,3 
follow the flows of the I r 's, and the action is the same as on the jm-torus, namely, 
S r = Aitj r . Along contour C4 the spinor rotation makes no contribution to the action 
while the half rotation around the torus in all three angles gives the action 



(82) 



As for the Maslov indices, we compute the complex matrix of Poisson brackets 
whose rows are indexed by the functions (I±, I 2 , I3, Jx, Jy, Jz), and whose columns 



are indexed by (zn,Zi 2 , z 2 i, z 22 , Z31, 
drop, the determinant is 
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(83) 



The final product of determinants is interesting, since these are the SU (2) invariants 
that Schwinger (Biedenharn and van Dam, 1965) used to construct the rotationally 
invariant state IJ1J2J3O) (with z rfl replaced by a\^). Bargmann (1962) and Roberts 
(1999) make use of the same invariants. For our purposes we need the winding number 
of the loop traced out in the complex plane by the product of the three determinants 
as we follow the four basis contours on the Wigner manifold. 

Proceeding as we did on the jm-torus, we find that the Maslov indices along the 
contours C r ,r — 1,2,3 are 4. For example, along the /i-flow z\\ and Z\ 2 get multiplied 
by exp(i , 0i/2), which causes two of the three determinants to be multiplied by the 
same factor, so the product gets multiplied by exp(ii/>i), which has winding number 
2 and hence Maslov index 4 when ipi goes from to Air. This is the same answer we 
found along the 7 r -flows on the jm-torus; this was not exactly a foregone conclusion, 
even though the contours are the same, because the tangent planes are different. On 
the other hand, in both cases we find the result (1801) for the eigenvalues ff u of the 
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operators 7 r , which of course must not depend on how we compute them. As for 
contour C4 the first leg, a rotation by angle 2tt about some axis, leaves all three 2x2 
determinants in (|83|) invariant, so the big determinant in the complex plane does not 
move. As for the second leg, since each ip r only goes from to 2ir we get a winding 
number of 1 along each J r flow, but since there are three of them the total winding 
number is 3 and Maslov index is 6. 

Combining this result with (|82| . we obtain the Bohr-Sommerfeld quantization 
condition for contour C4 in the form 

27r^j r — 6^ = 2n x integer, (84) 

r 

or, with J5DD, 

J2tf u = integer. (85) 

This is precisely the condition that the three quantum angular momenta must satisfy, 
in addition to the triangle inequalities, that they may add up to zero. It emerges in a 
semiclassical analysis because the Wigner manifold is not quantized otherwise. 

In conclusion, the Bohr-Sommerfeld quantization conditions applied to the jm- 
torus and the Wigner manifold give us a complete (and exact) accounting of all the 
quantum numbers and the restrictions on them that appear in the coupling of three 
angular momenta with a resultant of zero. It also allows us to identify the classical 
manifold (that is, its contour values) with a given set of quantum numbers. 

11. The amplitude determinant 

The generic semiclassical eigenfunction of a complete set of commuting observables is 
given by (TTJ) , with the amplitude determinant expressed in terms of Poisson brackets 
by ([2]). These formulas apply in particular to the state |iii2j3Wim2m3) on the left 
of the matrix element (|22|) . which is supported by the jro-torus in the large phase 
space. The state on the right, IJ1J2J3O), however, which is supported by the Wigner 
manifold, is an eigenfunction of observables that do not commute. Therefore we 
must rethink the derivation of Eqs. ([T|) and to see what changes in this case. 
In particular, we must see what happens to the Poisson bracket expression for the 
amplitude determinant, which is the solution of the simultaneous amplitude transport 
equations for the collection of observables. As it turns out, nothing changes, the wave 
function is still given by Eqs. (H|) and @, with the (now noncommuting) observables 
used in the amplitude determinant. In addition, there is a certain understanding about 
how the volume V in {1} is computed, since the angles a conjugate to the A's are no 
longer meaningful. 

Once this is done, we must evaluate the scalar product of the two wave functions 
by stationary phase. If both states were eigenstates of complete sets of commuting 
observables, then the answer would be ([5]) with amplitude determinant (j4|) , but again 
we must rethink the derivation of this result since the observables for one of the wave 
functions do not commute. Again, the answer turns out to be given by formulas (O 
and (0| of Sec. fJl with a proper understanding of the meanings of the volume factors. 

Having established these facts, we can then proceed to the (easy) calculation 
of the amplitude determinant for the 3j-symbol in terms of Poisson brackets, and 
finally put the remaining pieces together to get the leading asymptotic form of the 
3j-symbols. 
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11.1. Amplitude determinant for noncommuting observables 

We begin showing that Eqs. (Q]) and are valid for the state | Ji J2J3O) , with a 
proper definition of the volume factors. The classical functions defining the Wigner 
manifold are (7i, I 2 , I3, J x ,J y ,J z ). Let us refer to these collectively as Ak, k = 1, . . . , 6, 
let us write x l , i = 1, ... ,6 for the configuration space coordinates instead of the 
notation (xn, . . . ,£32) used above, and let us adopt the summation convention. The 
functions Ak form a Lie algebra, that is, {Ak, Ai} = cjj] A m , where c£J are the structure 
constants. The Wigner manifold is a compact group manifold with this Lie algebra, 
on which the Haar measure is both left- and right-invariant. This density is also 
invariant under the flows generated by the right-invariant vector fields, which in our 
case are the Hamiltonian flows of the functions Aj. The projection of this density 
onto configuration space is the density that provides the solution of the simultaneous 
amplitude transport equations for the functions Aj. These are the basic geometrical 
facts, which we now present more explicitly in coordinate language. 

The amplitude transport equations for the functions Ak, k = 1, . . . , 6 are 

— \pJx) 9Ak 
dx % dpi 



0, (86) 



where pi are the momenta conjugate to x l . These are six simultaneous equations 
that must be solved for the density SI (2;) on configuration space. Notice that 
dAk/dpi = {x l ,Ak} = the latter being notation we shall use for the velocity 

in configuration space along the Hamiltonian flow generated by Ak- The amplitude 
transport equation is a continuity equation, which is form-invariant under general 
coordinate transformations. 

Let us pick one of the branches of the inverse projection from configuration 
space onto the Lagrangian (Wigner) manifold. We shall suppress the branch index 
in the following. Let it 1 , i — 1,...,6 be an arbitrary set of local coordinates on 
the Wigner manifold, which we extend in a smooth but arbitrary manner into some 
small neighborhood of the Wigner manifold, so that partial derivatives of the u 1 with 
respect to all phase space coordinates are defined. Assuming we are not at a caustic, 
the transformation from x l to u l is locally one-to-one, and the Jacobian du l jdx^ is 
nonsingular. Under the inverse projection or coordinate transformation x — > u, the 
flow velocity transforms according to 

which defines the quantities XJ k \ . As a matrix, XJ k -> is nonsingular because the 
flow vectors are linearly independent on the Wigner manifold. As for the density, 
it transforms according to 

dx l 
du m 

so that the amplitude transport equations, lifted to the the Wigner manifold, become 
l-[o-(x)Xl k) } = 0. (89) 

Now define A^- as the matrix inverse to » 

^4)=^. (90) 



o~(u) — Sl(x) 



det • 
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As we will prove momentarily, the solution of Eqs. 

a{u) =|detAf ) |, 
which, by (|88|) , gives us the solution of ([86]) . 

. (fc) du l 



(91) 



det k i ( A 
detfej [ u 



dx 1 

dx 1 
{k) ~du~ l 



det fe z ( Xfa — 



dx 1 



detki x| fc ) 



= \det k i{x l ,A k }\ .(92) 



In carrying out these manipulations it is important to note that du l /dx° is taken at 
constant Ak, not pk- Thus the amplitude determinant for the wave function tp(x) 
associated with the Wigner manifold has the same Poisson bracket form shown in © , 
that is, in spite of the fact that the Ai do not commute. 

The essential differential geometry of these manipulations is that Xha = 
Xfyd/du 1 are the Hamiltonian vector fields on the Wigner manifold associated with 

functions Ak, A^ = A- fe ^ du l are the dual forms, X^X. 



(I) 



5 k , and a 



a(u) du l f\. . .Adu 6 is the Haar measure. The condition is equivalent to Lx (k) c = 0. 



To prove (|91|) in coordinates we substitute it into 
derivative, obtaining an expression proportional to 



and expand out the 



X) 



.-A 



(m) 



using commas for derivatives. Then we use the Lie bracket of the vector fields X^, 

— X\k)X(m)A ~ X (m) X (k),i = _c fem^(n)i (94) 



[X(k) , X( r 



where c£ m are the structure constants. Here we use the identity expressing the Lie 
bracket of Hamiltonian vector fields for two functions in terms of the Hamiltonian 
vector field of their Poisson bracket, [Xh,Xk] — —Xs HjK \ (Arnold, 1989). Thus 
(|93p becomes simply c™ m , which vanishes since for the group in question the structure 
constants are completely antisymmetric. 

Finally to normalize the semiclassical eigenfunction supported by the Wigner 
manifold we use the stationary phase approximation to compute the integral 

2 



dx 



br 



n(x) 1/2 exp[iS(x) - ifin/2] 



(95) 



where the sum is over branches and the branch index is suppressed. Cross terms do 
not contribute, and when the integral is lifted to the Wigner manifold it just gives the 
volume of that manifold with respect to the Haar measure, 

dx f 

= / du<j{u) = V w , (96) 



E 

bi 



|det{a;\ A,-}| 



where Vw is given by ([46 



11.2. Matrix elements for noncommuting observables 

Now we write the 3j matrix element (|22|) as (b\a), where A = (I\ , I2 , 13 , J x , Jy , Jz) and 
B = (Ii,l2,h, J\z, J2z, J-iz)- Actually this is not the most convenient form, since J z in 
the j4-list is a function of the J rz in the S-list. We fix this by performing a canonical 
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transformation [<f>x, </> 2 , <t>3, J\ z , J27., Jzz) ~^ (<Pi, 4>2, 03, J\z, J2z, Jz) on the functions in 
the -B-list, generated by 

^2(01, 4>2, 4>3, Jlz, J2z, Jz) = 4>lJ\z + §2-hz + 4>z(Jz ~ Jlz - Jiz)- (97) 

This gives J lz = Jlz, -hz = hz, Jz = Jiz + Jzz + Jsz, and <j> x = cj>i - cj> 3 , <f> 2 = (f>2 - 03, 
03 — 03- The linear transformation in the angles has unit determinant, so the 
volume of the jm-torus is still given by (f4*3"j) . Dropping the tildes, the i3-list is now 
(ii, J2, ^3, Jiz, J2z, Jz), which has four functions in common with the ^4-list. 
Now the integral we must evaluate is 

1 f 1 1 

(b\a) — —j^=^= I dx- 



VVXV^J Idet^,^}! 1 ^ \det{x\B,}\^ 



x 2^ exp{i[S A (x) - S B {x) - mtt/2]}, (98) 

br 

where the sum is over all branches of the projections of the two manifolds, and where 
fi just stands for whatever Maslov index appears in a given term (different /i's are 
not necessarily equal). An integral like this was evaluated by Littlejohn (1990), using 
the angles conjugate to the A's and B's, but those do not all exist in the present 
circumstances and we must evaluate the integral in a different way. 

Let us write A = (C,D) and B = (C,E), where C — I2, 13, J z ) are the 
four observables in common in the A- and i?-lists, and where D — (J x ,J y ) and 
E = (Ji z ,J2z) are the two pairs of observables that are distinct. The stationary 
phase set of the integral (|98|) consists of points x where 8(Sa — Sb)/9x 1 = 0, that 
is, it is the projection onto configuration space of the intersection of the A-manifold 
and the i?-manifold. That intersection, which we denote by /, was studied in Sec. [8] 
(it is a 4-torus). It is the simultaneous level set of all of the A's and B's, and at the 
same time the orbit of the commuting Hamiltonian flows generated by the C's. Its 
projection onto configuration space is a 4-dimensional region. 

We introduce a local coordinate transformation in configuration space x — > (y, z) 
where the four y's are coordinates along the stationary phase set and the two z's are 
transverse to it. We let the stationary phase set itself be specified by z = 0. We let 
(11, v) be the momenta conjugate to (y,z). Then the two amplitude determinants in 
(|98p may be combined with the Jacobian of the coordinate transformation to result 
in the square root of the product of two determinants, one of which is 

and the other of which is the same but with the substitutions A — > B, D — > E. But 
since the C's generate flows along /, we have {z\ Cj} — 0, and the lower left block of 
the two matrices vanishes. Thus, the product of the two determinants becomes 

[det{y, C}] 2 det{z, D} det{z, E}, (100) 

the square root of which appears in the denominator of the integrand. Evaluating the 
final two Poisson brackets in the (y, z; u, v) canonical coordinates, we have 

{z \ Dj } = ^i, {z i, Ej } = ^i. (101) 

We perform the z-integration by stationary phase, expanding Sa and Sb, 
regarded as functions of (y, z), to second order in z for a fixed value of y, and simply 
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evaluating the amplitude at z — (that is, on /). To within a phase, the z-integration 
gives 

-1/2 

(102) 



2n 



det 



d 2 S A 8 2 S B 



dzdz 



dzdz 



The determinant in this result must be multiplied by the determinants of the matrices 
(|101|) to get the overall determinant in the denominator after the z-integration. The 
product of these three determinants is the determinant of the matrix 

T 



dD 

dv 



dv\ 



dE 
dv 



(103) 



yCD \"'-/yCE\ 

where it is understood that a partial derivative stands for a matrix whose row index 
is given by the numerator and column index by the denominator, unless the matrix 
transpose or inverse is indicated, in which case the rule is reversed. Also, if a partial 
derivative is shown without subscripts, then it is assumed that it is computed in 
the canonical coordinates (y, z;u,v), and otherwise the variables to be held fixed 
are explicitly indicated. In the two middle matrices in (|103[) . the variables held 
fixed amount to differentiating v with respect to z along the A- and £>-manifolds, 
respectively, since dSA/dz = v(x,A) and dSs/dz — v(x,B). Notice that these two 
matrices are symmetric. 

Now we express the two matrices in the middle of (|103l) purely in terms of 
partial derivatives computed in the canonical (y, z] u, v) coordinates. We do this by 
writing out the Jacobian matrix d(y, z;C, D)/d(y, z;u,v) and the inverse Jacobian 
d{y, z; u, v)/d(y, z;C,D), multiplying the two together to obtain a series of identities 
connecting the forward and inverse Jacobian blocks, and then solving for the inverse 
Jacobian blocks in terms of the forward ones. We note that the block dC '/dv of the 
forward Jacobian vanishes, since it is {z, C}. Thus we find 



fdv^ 




\di, 


^ yCD ( 


/dv^ 




[dz, 





dD 

dv 



dC_ 

dz 



'dD | 


'dCY 


du ' 


K du) 


(T 


n 


V du , 





dC 

dz 



dD 

dz 



dE 
du 



dE 

dz 



dE 
dv 



-IT 



(104) 



Upon substituting these into (|103[) . that matrix becomes 

dD fdE\ T dD fdE\ T dD 
dv \dz J dz \dv J du 

T / rs^,\ -IT / ~ ^\ T 



dC 
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dC 
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dE 
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dC_ 

du 



dE 

du 



(105) 



where the first two terms are the beginning of the Poisson bracket {E, D}. As for the 
last two terms, we write out the vanishing Poisson brackets {C, E} and {C, D} in the 
(y, z; u, v) coordinates, making use of dC /dv — 0, to obtain 
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Actually the matrix of Poisson brackets {C, D} does not vanish everywhere in phase 
space, just on the A- (or Wigner) manifold, and in particular on the intersection / 
which is where we are evaluating them. Now substituting Eqs. (|106[) into the last two 
terms of (|105|) . those terms become 

(107) 
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in which the first two terms give us the remainder of the Poisson bracket {E,D}. 
As for the last major term, the factor in the square brackets vanishes, as we see by 
writing out the vanishing Poisson bracket {C, C} in coordinates (y, z; u, v) and using 
dC/dv = 0. 

As a result the integral (f9"5| becomes 



(b\a) = 



2tt 



VvXVb 
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br 
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|det{j/,C}||det{B,Z>}|V2' 



(109) 



where the branch sum runs over all branches of the projection of I onto configuration 
space as well as the two disconnected components of I (the two 4-tori discussed in 
Sec. [8]), and where Si is the phase on a given connected component of / (this is the 
phase ±Sj m computed in Sec.[9|). We have also dropped an overall phase, and we are 
not attempting to compute the Maslov indices in detail. The amplitude determinant 
has been reduced to a 2 x 2 matrix of Poisson brackets of the observables in the A- 
and B- lists that differ, exactly as in ([4}. Calculating this matrix explicitly, we find 



\&et{E,D}\ 



{Jlzi Jx\ {Jlz->Jy\ Jyl Jxl 

{J2z,Jx} {J2z,Jy} Jy2 ~Jx2 

= \JxlJy2 ~ Jx2Jyl\ = |z • (J X X J 2 )| = 2A Z , (110) 

where A z is the projection of the area of the triangle A onto the x-y plane (see Eq. ([50| 
and Fig. 2 of Ponzano and Regge (1968)). This quantity is invariant under rotations 
about the z-axis, that is, it Poisson commutes with J z . It also Poisson commutes 
with the other three variables in the C-list, (Ji, 12,1s), and so is constant on the 
intersection / and can be taken out of the y-integral. The same applies to the phase 
factor, since Si is also constant on the /-manifold. Then the y-integral can be done, 
since | det{y, C}\ is just the Jacobian connecting y with the angle variables conjugate 
to C — (ii, I2, h, J z ), denoted above by (ipi,4>2, 4>3i4')- Thus the y- integral just gives 
the volume Vi of the intersection / with respect to these angles, see In fact, 

had the variables C not been commuting, but if they had formed a Lie algebra, then 
Vi would be the volume of / with respect to the Haar measure of the corresponding 
group. This circumstance arises, for example, in a similar treatment of the 6j-symbol. 

As a result of these rather lengthy manipulations of amplitude determinants, we 
obtain the final, simple result, 

2tt ^-^ gi(Sj-/i7r/2) 
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where the branches now run over just the two disconnected pieces of the intersection 
/. This is a version of (|5|), with the right understanding of the volume measures, 
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generalized to the case at hand in which the observables do not commute. The actual 
calculation of the final amplitude determinant takes just one line, Eq. (|110[) . 

In fact, for our application the volume Vj and the remaining amplitude 
determinant are the same for both branches and can be taken out of the sum. The 
relative Maslov index between the two branches is 1; we will not belabor this point 
since the answer is already known. We simply note that by splitting the Maslov phase 
in/2 between the two branches and subsitituting Va — Vw, Vb = Vj m , we obtain to 
within an overall phase the result of Ponzano and Regge, 



12. Conclusions 

In many ways the 3j-symbol is not as interesting as the 6j-symbol, of which it is 
a limiting case. We intended our work on the 3j-symbol as a warm-up exercise, 
expecting a routine application of semiclassical methods for integrable systems. 
The nongeneric Lagrangian (Wigner) manifold was a surprise. Similar nongeneric 
Lagrangian manifolds occur also in the semiclassical analysis of the 6j- and 9j-symbols. 

If all one wants is a derivation of an asymptotic formula, then there are many 
ways to proceed. For example, one can simply take the expression for the symbol 
due to Wigner (3j) or Racah (6j) as a sum over a single index, and apply standard 
asymptotic methods (Stirling's approximation, Poisson sum rule, etc). But if one 
wants a derivation that reveals the geometrical meaning of the classical objects that 
emerge (the triangle, the tetrahedron, etc), then an approach such as ours may be 
preferable. 

Our approach is more geometrical than earlier ones, and in that respect is closer in 
spirit to the work of Roberts (1999), Freidel and Louapre (2003) and later authors. It 
is likely that at some deeper level all these methods are the same, although superficially 
we see only a little similarity between our work and these others. 

One may also desire a method that makes the symmetries of the symbol manifest. 
Our analysis does not do this for the 3j-symbol, but those symmetries are not manifest 
in Wigner's definition of the 3j-symbol that we employ as our starting point, either. 
To bring the symmetries out it seems necessary to employ some construction related 
to Schwinger's generating functions, which involve lifting the definitions into higher 
dimensional spaces. 

Our method of calculating amplitude determinants in terms of Poisson brackets 
may have computational advantages in other applications, as well. The method can 
be remarkably easy to use. For example, the 6j-symbol can be defined as a matrix 
element, 



which is the unitary matrix in (j'12, J23) defining a change of basis in the subspace in 
which four angular momenta of given lengths add up to zero (0 means J = 0). In this 
case there are eight observables on each side of the matrix element, of which seven 
are common and one is different. Thus the amplitude of the 6j-symbol is the inverse 
square root of the single Poisson bracket, 




(112) 




(113) 



{j| 3 ,J? 2 } = 4J 1 -(J 2 x J 3 ), 



(114) 
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as follows immediately from (|37p . One sees immediately that it is proportional to the 
volume of the tetrahedron. A similarly easy calculation is possible for the 9j-symbol. 
It is harder, however, to express these amplitudes in terms of the quantum numbers 
(the magnitudes j r ), that is, to translate these magnitudes into vectors 3 r that lie on 
the stationary phase set. We shall report on these and other extensions of our work 
in future publications. 
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Figure 1. The action of an SU(2) rotation about a fixed axis on a point of the 
large phase space, and its projection onto angular momentum space. 



Figure 2. Harmonic oscillator motion in the x rll -p rlJ , plane, generated by I rfl . 
Definition of angle 8 rll is shown. 



Figure 3. The quantum state \j\32jzm\ m2m^) corresponds classically to a 3- 
torus in the small phase space S 2 X S 2 X S 2 , which can be visualized as three J 
vectors lying on three cones with fixed values of J z . The three azimuthal angles 
are independent. 



Figure 4. If ji, j% and j'3 satisfy the triangle inequalities, then they define a 
triangle that is unique apart from its orientation. A standard orientation places 
the triangle in the x-z plane with sides Ji, J2, J3 oriented as shown. The angle 
opposite J r is r\ T . 



Figure 5. An schematic illustration showing how the Wigner manifold in the 
large phase space is the inverse image under 1 of a set of triangles formed 
from three angular momentum vectors with vanishing sum, all related by rigid 
rotations. 



Figure 6. By rotating the reference orientation of the triangle about the y-axis, 
we can give J3 the desired projection 1113 onto the z-axis. 



Figure 7. Once vector J3 has the desired projection ms, we rotate the triangle by 
angle 7 about the axis J3 to make J 2 have its desired projection m^. This cannot 
always be done for real angles 7, but when it can be done there are generically 
two angles that work, illustated by points Q and Q' in the figure. 
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